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I. INTRODUCTION 


a matrix of the 

form 


/ Xq Xn-I 

• ■ • X2 

Xl 

Xl Xq 

■■■ X3 

X2 

X2 Xl 

■ ■ ■ X4 

X3 

XN-2 Xn -3 

■■■ Xq 

Xn-1 

\ Xn-1 Xn-2 

■■■ Xl 

Xq 

as are circular 

permutations 


= CIRC [xq,xi, ... ,xn-i\ 


authors is to take successive rows as circular permutations of the first row.) It is a particular form of a Toeplitz 
matrix, a matrix whose elements are constant along the diagonals. The N eigenvalues of the N x N circulant matrix 
above are 


Xq + UJ^Xi + + ■ • ■ + w 




XN-1, p = 0,...,iV-l 


( 1 ) 


where w = is an A^th root of unity. The determinant of the matrix is therefore the product of these eigenvalues: 


Xq 

Xn-1 ■ 

■ X2 

Xl 

Xl 

Xq ■ 

■ X3 

X2 

X2 

Xl ■ 

■ X4 

X3 

Xn-2 

XN-3 • 

■ Xq 

Xn-1 

Xn-1 

Xn-2 ■ 

■ Xl 

Xq 


N-l 


= det[xo,a:i,... 


Xq + OJ^Xi + • • • + W 


(JV-l)p 


XN-1 


P—0 




Some examples of these determinants for small values of N are: 
det)^,^,^] = + -?,ABC-, 


( 2 ) 


det[A, B, C, D] = - B'^ + - 2A^C‘^ + 2B‘^D‘^ - AA^BD + AAB'^C - ABC^D + AACD^; 


det[A,B,C,D,E] = A^ + B^ + 

-bA^BE - bA^CD - bAB^C - bB^DE - bAC^E 
-bBC^D - bABD^ - bCD^E - bADE^ - bBCE^ 

+bA^B^D + bA^BC^ + bA^CE"^ + bA^D'^E + bAB^E^ 

+bAC^D^ + bB^C^E + bB'^CD'^ + bBE^E"^ + bC'^DE'^ 

-bABCDE. 

The product on the right in eq.(2), when expanded out, contains terms. There are considerably fewer terms 
in these examples, and show some apparent patterns in the coefficients. There appears then to be considerable 
cancelations occurring in the expansion of eq. (2), as there must be, since all the complex terms cancel. This suggests 
that there is a simpler way to calculate these determinants. In addition, expression (2) has the disadvantage of using 
a complex number, w, to find integer coefficients. 


This work was initiated then by the desire to derive an explicit formula for the coefficients in the expansion of the 
determinant, and to see if the apparent patterns in the coefficients for small N appear in the general case as well. 

Near the completion of these results I discovered the monograph on circulant matrices by Wyn-jones [l|, where 
chapters 10 and 11 are devoted to the same problem that inspired this effort, that of finding a formula for the 
coefficients. The formula in ref. [1], which was first stated in a 1951 article by Ore [ 3 ], bears some similarities to the 
one derived in this article; both formulas involve sums over the partitions of sets of integers, for example. But there 
are significant differences between the two, and Wyn-jones’ proof uses a completely different method than the one 
presented here. In addition, it is likely that one formula is more efficient than the other for some coefficients and less 
efficient for others. So it seems reasonable to present this alternative one as well. Some of the other results in this 
article were proved in ref. [1], but again the proofs here are by different methods. 
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II. PRELIMINARY RESULTS: SYMMETRIES OE THE Ur^, COEFFICIENTS 


As shown in the examples in Section I, the determinants can be expanded as sums over products of the matrix 
elements: 


det [a:o,a;i,... ,xn-i] = E ^“0 ■■■aAT-i^ao 

O^dO — — 1 

The sum can also be written in the form 

det [xo,xi,.. .,XN-i] = E 

0^ A^o 5 ■ ■ — 1 


Mo 


Mn-1 
^N-1 ’ 


(3a) 


(3b) 


where C’mo -Mjv-i “ Cag-.-aN-n ^0 + ••• + Mff-i = N, and Mm is the multiplicity of m, the number of times 
the integer m occurs in the index set [oq, ... ,a 7 v-i]- Let >5[ao. ',aiv-i] permutations of the index set 

[oo,..., uat-i] : cr[ao,..., oat-i] —>■ [gq, ..., on-i]- From eq.(2), the coefficients in eq.(3a) are 


Ca 


E 


OJ 


(7i+2(T2H- \-{N — 1)(7N-1 


( 4 ) 


where the sum is over all elements of 5[ao,... ,ajv-i] ■ One could of course use this expression to calculate these coeffi¬ 
cients, but it is not very convenient, involving as it does a sum over a large number of complex numbers. 

We will use the notations [a] = [ao)'"' jOat-i], (oq < Qi < • • • < oln-i)^ and [a] = [cti, ..., ctat-i], where cr is an 
arbitrary permutation of [a]. Let be the subset of <S[o] defined by 

^[a] = { cr € 5[a] I CTi -b 2 ct 2 H-h (A - l)crAr-i = k mod N } 

A(k) (k) (k) 

and let be the number of elements in , which we will denote as . More generally, if 5 is a finite set, 

then is the number of elements in S. The elements of OjO will be referred to as “fc-mod permutations” (of [a]). 
Since ^ eq.(4) reduces to the form 






with 




■A' 


(N-l) 


m 


Affp! • • • Ad AT—1! 


( 5 ) 


( 6 ) 


When all the a’s are equal, (oq = oi = • • • = oat-i = a; Ma = N), 4S[aa- -o] Las just one element, so there is only one 
term in the sum in eq.(4). The corresponding coefhcient is 

_ 27rio(l+2H \-N — l)/N 7 ria(N—l) / -|\a(Ar—1) 

^o---a — C — e — \k) 

In the rest of this article we will consider coefficients with more than one multiplicity. 

We have = 1; more generally, (w"^)^ = 1 for any integer d. If d ^ 0 mod N, the set ..., = 

1} consists of the zeros of the polynomial — 1: 

Z^-1 = {z-UJ^){z-UJ^^)---{z-UJ^^-^^^){z-l). 

Expanding on the right, we have 

and since the coefficient of on the left is zero we get the identity, 

N-l 

9=1 

where the mod-Y Kronecker delta is defined as 
1 if d = 0 mod N, 


( 7 ) 


zN _ 

^f/,0 — 


0 otherwise. 
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Theorem 1 

C[a] =0 if Go + ai + ■ ■ • + gn-i ^ 0 mod N. 


(This is Proposition (10.4.3) in ref.[l].) 

Proof: Let [o] be such that 

ao + ai + • • • + oat-i = X mod IV, X ^ 0, 
and therefore, for any permutation [a] of [a], 

(To + ai + ■ • ■ + gn-i = X mod N. 

We can assume X to be positive; (if not, we can make the replacement X ^ X' = N — \X\). If we add the second 
equation repeatedly to the condition 

(Ti + 2 ct2 + ■ ■ ■ + {N — l)(TAr-i = k mod N, 

for permutations in we get 

(To + 2 (Ti + ■ • ■ + {N — l)aN -2 = {k + X) mod N, 
gn —1 + 2(To + • • • + {N — l)(TAf —3 = (fc + 2X') mod N, 


(72 + 2(73 + ■ ■ ■ + {N — l)(7l 


{k + N — X) mod N. 


There is therefore, for each k, a sequence of one-to-one, onto mappings: 

[a] [a] [a] [o] [a] 

induced by circular permutations oi gq, ■ ■ ■ ,gn-i- Accordingly, 

/ll^) _ A(f^+X) _ Ak+2X) _ Ak+N-X) r 7 _ Q 1 X — ^ 

^fal — ^[o,l — ^[o,l — — ^[o,l J lor K — U, i, . . . , ^ i. 


If A = 1, or if A and N are relatively prime, then = .. • = and the summation in eq.(5) for 

C[a] equals zero as a result of the identity (7). If on the other hand gcd(A, A) > I, then MX = 0 mod N for some 
M < N, and we have instead A sets of equalities: 


A 

A 


( 0 ) 

[a] 

( 1 ) 



= A 


(i+x) 


.((M-l)X) 

^[a] 

_ .(1+(M-1)X) 


7(X-1) _ 4(2X-1) _ _ .(X-l) 


In this case, eq.(5) becomes 


q.] = ( +-1(1’“+- ■+)(1++ ■ -+). 

But since oj^ = fQj- some integer x ^ 0 mod N, the sum 1 -|- oj^ -|- -I- • • • -|- equals 

zero as a result of the identity corresponding to eq. (7) with the substitutions N ^ M and w —>■ wm = . A 


A number of corollaries follow from Theorem 1. 

Corollary 1 For even N, Coi23---n-i = 0. 

Corollary 2 For odd N, Cao---aN-i = 0 */ the multiplicities are 2, 1, ..., 1, 0. 
Corollary 3 Ca-.-ab -b = 0 if the multiplicities Ma and Mh are relatively prime to N. 
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The proofs are straightforward and we omit them. 

As noted, if [oq, ..., oat-i] satisfies the condition, 

Oq + Oi + • • • + otv-i = 0 mod N, (8a) 

then any permutation [a] of [a] also satisfies it. Expressed in terms of the set of multiplicities, [M] = [Mq, Mi, ..., 

0 < Mo, Ml,..., Mm -1 < N, this condition becomes 

0 • Mo + 1 • Ml + • • • + (A^ - 1) ■ Mn -1 = 0 mod N, (8b) 

subject to the restriction 

Mq + Ml + • • • + Mat—1 = N- (8c) 


Since we are only interested in non-zero coefficients, unless otherwise stated we will assume that all [a]’s and [trj’s 
satisfy (8a), and all [M]’s satisfy (8b,c). 

We define P to be the circular permutation operator, P((Jq) = <Jq+i'. 

P[ct] = P[cto, (Ti, . . . , crAr-2, O-AT-l] = [cti, CT 2 ,..., crAr-2, O'AT-l, CTo]. (9) 

Let {a} = {(7o,..., ctat-i} denote the equivalence class under P of a /c-mod permutation [a\. 

{ct} = distinct elements of the set { [cr], P[cr],..., P^~^[a] }. 

If [a], with multiplicities Mq, Mi, ..., Mjv-i, is such that gcd(Mo, Mi..., Mat-i) = 1, each equivalence class of the 
permutations of [a] has N elements and each is N times the total number of such equivalence classes. Therefore 

(k) 

C\^a] = 0 mod N. On the other hand, if the multiplicities have a common divisor d and if a permutation in 
takes the form [upaqUr • ■ • Ug Op • ■ • Ug apUgUr ... a^] where the sequence apUgUr • ■ • Ug is repeated d times, then the 
equivalence class containing this permutation has only N/d distinct elements. So, from eq.(5), we have 

Proposition 1 If gcd{Mo, Mi,..., = d, then Cao---aN-i =0 (^/d)- 

For an integer n and a permutation [ uq, ..., on-i ] 'we define the operations 

n+ [ (To, ■ • ■ ,0'jv-i ] = [ n-I-(To,... ,n-I-CTjv-i ], 

n X [ (To,..., (Tat-i ] = [ UCTo,..., ncTAT-i ], 

where addition and multiplication is modulo N. Further, we define |(t| to be the permutation [o'] of [(To, ai,..., ctat-i] 
that puts it into ascending order: (Tq < < • • • < 

Proposition 2 

(A) Caq...a^_, = (-l)”('^-l)q„+[ao.....a„_q|. 

(B) Cao-aN-i = C'|nx[ao-ajv-i]| M ^ o-ud N relatively prime.. 

(This is Proposition 10.5.5 in ref [1].) 

Proof: 

(A) Condition (8a) is invariant under the transformation Uk {n + au) mod A^ for n = l,...,Af — 1, while 
CTi -I- 2 ct2 {N — l)CTAr-i changes by 

A (cTi -I- 2(T2 -I- V{N- l)(TAr_i) = n ' 

This transformation will then give a factor of to each term on the right side of eq.(4), so 

/~i _ ( _-|\n(Ar—1)^ 

'-^|n+ao,---,n+aAr-i I V '-^ao---aN-i • 
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(B) The set contains the same elements, up to order, as the set {l,a;,a;^,... when 

N and n are coprime. We can therefore make the replacement w uj'^ in all the terms and factors in the product on 
the right in eq.(2), since this amounts to just a re-ordering of the factors. Eq.(4) then becomes 


Q, 


E 




+2(72 H- \-{N—1)(t n — 1 


(TGS\r 


E 


nai+2nrT2-\ -h(A''-l)rao-jv-i _ , ,, 

^ — (-'|(nx [a] I ■ 


A 


In section IV below, we will use the two symmetry operations above to classify the coefficients, (or, equivalently, 
their multiplicative index sets), into “additive multiplets”, in which the coefficients are related by the additive sym¬ 
metry, and into “super-multiplets”, in which they are related by the additive symmetry, the multiplicative symmetry, 
or a combination of both. We will use the notation {C[a]} or to represent the additive multiplet containing 

(7[a], and (Cfaj) or to represent the corresponding super-multiplet. As a result of Proposition 2, coefficients 

in an additive multiplet or in a super-multiplet are equal, up to sign, and have the same multiplicities, up to order. 
The number of super-multiplets is thus equal to the number of coefficients that need to be calculated to solve for the 
determinant. 


As pointed out by Wyn-jones, condition ( 8 a) is a necessary but not sufficient condition for a coefficient to be 
non-zero; he gives two examples, (^001335 and CooDOiiiseSj which are zero despite satisfying ( 8 a). A sufficient (but not 
necessary) criteria for “condition-( 8 )” coefficients to be zero is given in Section HID, as a corollary to Theorem 3. 


Wyn-jones did prove that, for N a prime number, condition ( 8 a) is a sufficient condition for C^a] to be non-zero. 
An alternate proof of this is given as a corollary of Theorem 2 below. 

The Mobius function Q ^(n) is defined as 

{ 1 if n = 1 ; 

0 if n has one or more repeated prime factors; 

(— 1 )^ if n is the product of k distinct prime factors. 

Then 

Theorem 2 



where, in the sum over the divisors of N, is to be identified with . 

Proof: We define the mapping n ■ [a] ^ [a'], where n is an integer that is coprime to N , by cfq ^ a'q = mod n, 
(or just (Tnq, with mod N multiplication understood). Under this map. 


Af-l N-1 N-1 N-1 

E E ^ E E 

(J=0 q=0 q'=0 (j'=0 

where n~^ is the multiplicative inverse of n. [a'\ is therefore an element of and n induces a mapping: 


n 


(k) 




{n ^k) 


If n is coprime to N then so is n which induces the inverse map: 


n 


(k) 


n 


(n ^k) 


This is therefore a one-to-one onto mapping between and These two sets have the same number of 

elements then, and Aj^j^ = Aj"*^ for all n coprime to N . In particular, Aj^j^ = Aj"j^ for all n coprime to N and, if N 
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is a prime number, = • • • = 

n ! [a] [a] [a] 

Let {do, di,... jdx} be the set of divisors of TV: I = do < di < ■ ■ ■ < dx = N. With the replacement ^ 
we separate the sum in eq.(5) in the manner below: 

7V-1 N-di N-dK-i 

gcd{N,q) = l gcd{N,q)=di gcd{N ,q)=dK-l 

In the sum corresponding to gcd(A^, q) = dp, the summation variable q equals kdp for some k coprime to N. By the 
relation above, for this sum. The A coefficients can therefore be brought outside: 

K-l N/dp-l K N/dp-l 

=-4™ + E <]■’ E E 

p—0 gcd{N,k) — l p—0 gcd(N,k) — l 

The result then follows by noting that the sum over k 


N/dp-l 
gcd(A^,fc) —1 

is the Mobius function /r(iV/dp)0. A 


The coefficient C^a] can thus be expressed as a sum of integers, without the use of complex numbers. This expression 

(k) 

is not very useful in calculating C[a] since it still leaves us with the need to count the number of elements in the 
sets, although for fewer of these sets than appear in eq.(5). 


Using this theorem though we can prove that C[a\ is non-zero when N is prime, provided of course that the index 
set [a] satisfies condition (8a). Changing back to the previous notation , we have 

Corollary 4 For N prime, 


Cr„i = 


N 


- 


(iV-1)! 


7V-1 V Mo\---Mm-i\ 


^0. 


Proof: The divisors of N consist of just 1 and N, so from Theorem 2 


Cr„i = - a: 


l(l) 


We also have in this case that 


N\ 


Mq\ ■ ■ ■ Mn—i\ 


= A 


(0) + + • • • + AYr^’ = +iN- 

[a] ' [a] ' [a] [a] '' '' [a] 


t(l) 




(0) 


l(l) 


from the proof of Theorem 2 and from eq.(6). The first equality follows from these relations. To show that C[a] ^ 0, 
Ca-a = ^ 0, SO we need consider just the cases where all the multiplicities are less than TV. For these, 

oc A, so = TVUjq] for some integer If C[a] = 0, then 

(iV-I)! 

N must therefore divide {N — I)!. But this is a contradiction since N is prime and all the factors in {N — 1)! are less 
than TV. A. 


Corollary 5 in Section IIIB below, also a corollary to Theorem 3, extends the criteria for C[a] ^ 0, i.e., TV being a 
prime number, to TV not dividing • 







III. DERIVATION OF A FORMULA FOR G 




A. Definitions and statement of the theorem 


We will denote the set of partitions of the integer phy Vp, with a particular partition consisting of j parts denoted 
by Z = {zi,..., Zj) or, equivalently, by the set of multiplicities [ki,... ,kp]: 

P = Zi-\ - \- Zj 

= ki 2/l 2 “t“ ■' ■ “t“ pkp] k\ ‘ kp = j. 

The set of partitions of a set of not-necessarily-distinct integers {a, 6 ,..., c} will be denoted by V{ah ■ ■ c}. For 
a particular partition 0 S V{ah • • • c}, 0 = ( 6 *i) • • • (0j), where each 0^ is a subset of {a, 6 ,..., c}. We will express 
these subsets, for example, as 6 * = abc, in contexts where there should be no confusion. As an illustration, the set of 
partitions of the set {a, a, &, c}, where a ^ ^ c, is 


7^{aa&c} = I (aabc), {aab){c), {aac){b), (a 6 c)(a), {aa){bc), {ab){ac), 

{aa){b){c), (a&)(a)(c), (ac)(a)( 6 ), ( 6 c)(a)(a), {a){a){b){c) | . 


For a partition 0 and for any subset 0 of {a, c}. Kg will denote the multiplicity of 9 in the partition; for 9 G Q, 

(n\ 

rria ' will denote the multiplicity of the element a in 9. Finally, the sum of the elements of 9 will be denoted as the 
trace of 9, Tr(0): 


Tr(0) = ^ a. 
aee 


The map given by Zk = sends the set partition 0 = (0i) • • • {9j) to the integer partition Z = {zi,. .., Zj). Two 
partitions 0 and 0' in V{ab- ■ -c} are equivalent iff they map to the same partition Z. In this case, each 9k and 
have the same number, zr, of elements. The 0 equivalence classes can then be labeled by the partition Z, and each 
0 e {Z} for some Z. 

We will also use the multinomial coefficient (p; fci,..., kp)* , defined as 0 


(p; ki,k2,...,kp)* 


_p!_ 

1 * 1 /ci! 2*2^2! • • -p^pfcp! ’ 


p = fci + 2^2 H- \-pkp, 


and the Heaviside step function 
( 1 if n > 0 , 

H[n] = 

I 0 if n < 0 . 


For convenience, we will write the coefficient Caa-.-an-i C'd-.-o i - i Ai---ApAp+i with the O’s and I’s separated off, 
with Ai = omo+Mi, Ap+i = ajv-i, and p = N — Mq — Mi — 1 > 0. The A indices lie between 2 and N — 1, with 
multiplicities M 2 , M 3 ,..., M^v-i- 

With these definitions. Theorem 3 is the statement below. 


Theorem 3 

The coefficients with multiplicities 0 < Mg, Mi,... M^v-i < — 1 in the expansion (3a) are given by 



9 


^ ( -\\N — Mq — 1-\t 

<-'0---0 1---lAi---ApAp+i — 0q,Mi+2M2 + --- + (N-1)Mn-i v 

{N -Mo- 1)! 1 


i: n^n 

0e7’{Ai-Ap} 9 


{zk - 1)! 


^ ^ ^ ^P + 1 OitzTJS A. ... A \ ft ^ lc=^ ^9 -! 


N-V 


X El-")' E 

M=1 Ai,...,Aj=0,l 


^MjAiH -hAj H 


Mi-^XtXiOt) 


t=i 


N — Mo — 1 — X/t=i ^t{X{9t) + zt) \ ^ / X{9s) + Zs — 1 


Ml - ELi ^tX{9t) 


n 


Zs -1 


where j is the number of parts in the partition 0, Zs = ffhg is the number of elements in the part 9s, and 
X{9s) = —Tr{9s) mod N 
such that 0 < X{9s) < N — 1. 


For each partition, the sum over the A’s has 2^—1 terms. However, the presence of the step function in the 
expression means that, for a term to be non-zero, it must satisfy the condition 


Mi-^XtX{9t) > 0. 


t=i 


As a consequence, for many coefficients, many if not most of the terms in the sum over the A’s are zero. 

As an example, and as a comparison to the method used in [l[, we will calculate the coefficient Cooiiii 3788 using 
this theorem. In this case, the k, m and the multiplicities all equal 1 except for Ma^^-,- We have then 


Cooiiii3788 — —10 ■ 


7! 


4!1!1!2! 




^M,AiH -hAj 


0e7’{378} 


M=1 


Ai.....Aj=0,l 


X H 


4 - ^ XtX{9t) 


7 — Et=i ^t{X{9t) + Zt) \ 1 / X{9s) + Zs — 1 


4 - ELi ^tX(9t) 


n 


Zs - 1 


. 


The partitions of the set {3, 7, 8} and the values of (zi ,..., Zj) and (Ai,..., Xj) for the non-zero terms are: 
(378) : zi=3; A(378) = 2; Ai = 1; 


(37) (8) : 

(38) (7) : 
(3)(7)(8) : 

Then 

Cooilll3788 


(zi,Z 2 ) = (2,l); A(37) = 0, A(8) = 2; (Ai, A 2 ) = (1,0), (0,1), (1,1); 

(zi,Z2) = ( 2 , 1 ); A(38) = 9, A(7) = 3; (Ai, A 2 ) = ( 0 , 1 ); 

(^i,^ 2 ,Z 3 ) = (1,1,1); W(3) = 7, A(7)=3, X(8)=2; (Ai, A 2 , A 3 ) = (0,1,0), (0,0,1). 



= -10{ 105 - 175 + 50 } = 200, 
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in agreement with the value obtained in ref. [1]. 

By way of comparison, the method used by Wyn-jones involves partitions of the partial index set {1,1,1,1, 3, 7, 8,8} 
that are make up entirely of “null multisets”; i.e. (in my notation) subsets 0 of {1,1,1,1,3, 7,8, 8} such that Tr(0) = 0 
mod N. Such a partition consisting of j subsets contributes to the term in the expansion, unlike the calculation 
above, where each partition can contribute to several terms. Although P(11113788) contains many more partitions 
than P(378), (4,140 compared to 5), the calculation of this coefficient in ref. [1] involves only the 5 partitions, 

(37)(118)(118), (1117)(1388), (37)(111188), (118)(11378), (11113788). 


B. Proof of Theorem 3 

We will assume in the following that the indices of all coefficients considered satisfy conditions (8) and will omit 
writing out explicitly the mod-A^ Kronecker delta _ ^ in the formulas for the coefficient. 

Consider first the case Ma = 1 and the coefficient Cq-.-o i-.-iyi- In this case we have Mq + Mi + 1 = and 
0 • Mo + 1 ■ Ml + A • 1 = 0, so A = Af - Ml. 

With the freedom we have from the circular-permutation symmetry we can choose, when counting the number of 
elements in each set fIo^.^.oi...i(iv_Mi)’ 1° consider only permutations of [0 ■ • ■ 0 1 • • • 1 (A^ — Mi)] for which 

ao = A = N — Ml, 

CTpi = Crp2 = • • • = =1, 1 < Pi < P2 < • • • < PMi < N -1, 

all other aq = 0. 

(k) 

The number of elements in i...i ^ times the number of such permutations that satisfy the condition 

k = 0 ■ ao + I ■ (Ji -\ - h (A^ - 1) ■ (JAT-i = Pi -I- p2 H - f PMi ■ 


^lyloi...i(N-Mi) therefore equals N times the number of partitions (mod N) of k into exactly Mi distinct parts less 
than or equal to A^ — 1. We define Q{n] X, b) to be the restricted, mod-A^ partition function 


Q{n-,X,b) = #{ (pi,p 2 ,... ,px) I Pi H-hpx = n mod N; l<pi<---<px<b< A^}; 

i.e., the number of partitions mod N oi n into X, distinct parts, with each part less than or equal to b. For 
n = k, X = Ml, and b = A^ — 1, each such partition corresponds to a circular-permutation equivalence class 
{N - ... ,aiv_i}. Co...oi-i(iV-Mi) is then, from eq.(5) 


N-l N-l 

C'o---oi---i(Ar-Mi) = X] ■ X Q{k;Mi,N - 

k^O k^O 

By inspection, it can be seen that the partition function Q{n; N — 1) has the generating function 

N-lN-l N-l 

^ ^ g(n;A,Af-l)p^w" = J] (1 + pw«) = 1+ ••• +y^(wPi+-+P»+•••) + ■•• . 

X—0 n—0 q—1 

Now referring back to eq.(2), the product over q in the equation above is, except for the missing q = 0 factor, the 
determinant of an A^ x N circulant matrix with elements 1, fc, 0,..., 0. Using this fact and evaluating the determinant, 
we find 


N-l 

P (1 -k yw«) 

9=1 


1 

1 + y 


10 0 ■••Op 

pi 0 ■••00 

Op 1 ■••00 


0 0 p 1 0 

0 0 0 ■••pi 


N-l 

= II (-y) 

m—O 


NxN 
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Then 

N-l 

^g(n;X,iV-lV" = (-l)^. 

n—0 

and our result for this coefficient is 

Co...01---l(Ar-Mi) = ( —1)'^^-/V. 

Next we consider the coefficient C'c-.-oi---! AiAa - ApAp+i, with p = N — Mq — Mi — 1. For given sets {pi,... ,PMi} 
and {qi,... ,qp} of, respectively, Mi and p distinct integers between 1 and — 1, we assign the permutation sending 
[0 ■ • • 0 1 • • • 1 Ai ■ • • Apj^i] to [(To • • • (Tat-i] as follows: 

(To = Apj^i] 

<^Pi= ■■■ = =1. 1 < Pi < • • • < PMi < - 1; 

(Tgi = Ai, Uq^ = A 2 , ag^ = Ap; 1 < qi, 92 , • ■ ■, qp <-N - 1, qi ^ q, for all i j. 

(k) 

The requirement that a particular permutation be an element of ^low 

Pi H-h PMi = k- Aiqi - A2q2 - Apqp. 

For the above set of distinct integers {qj}, 1 < qi, ■ ■ ■ ,qp < N —l,we define Q{n-, X, 6 | < 71 ,..., (jfp) to be the mod-A^ 
restricted partition function. 


Q{n;X,b \ qi,...,qp) 

= # I (Pi, ■ • ■ ,Px) I Pi H- \-Px=n mod N; I < Pi < ■ ■ ■ < px < b; Pi ^ qj for all i and j |. 

As before, Q is the number of partitions mod N of n into X, distinct parts less than or equal to b, but now with the 

(k) 

additional restrictions that none of the parts are equal to any qj. The set of all permutations in flg.;.g ^ is 

then the sum of N ■ Q{k — Aiqi — ■■ ■ — Apqp] Mi,N — 1 | qi,... ,qp) over all allowed values of the g’s, divided by the 
factorials of the multiplicities to avoid overcounting: 


"^0---0 1---lAi-"Ap+i 


N 

ibg! • • • Mx-i^- 


Q{k-Aiqi - 

l<qi,...,qp<N-l 


Apqp] Ml, N 


1 I qi, ■■■,%) 


qi¥^qi 


Le. the sum is over all lattice points {qi,... ,qp), ( 7 i € [1, A^ — 1] such that no coordinates are equal. In the following, 
we will omit explicitly writing out the limits of the summation over the q variables, with the understanding that the 
limits will, unless otherwise indicated, always be from 1 to A^ — 1. 


[Note: As pointed out in a previous section, in the case where the multiplicities of an index set all have a com¬ 
mon factor d, there are some equivalence classes that have only N/d distinct elements. The expression above is 
however valid in those cases as well; see Appendix B for a discussion of this issue. ] 

Therefore, again using eq.(5). 


C'o...O l -l Ai -Ap+i = 


N 


N-l 


M21---Mn-iI 


N 


Q{k - Aiqi - Apqp] Ml, N - l\qi,...,qp) 


k—0 Qi^qj 


N-l 


Af^2' ■ ' ■ ^N—1 • 


^A^q,+-+Apqp Y Q{k']Ml,N-\\qi,...,qp) 


qi¥^qi 


fc '=0 


(10a) 


where we’ve changed the summation index to k' = k — Aqi — ... — AqMA- 1'1 since Q uses modular arithmetic, this 
does not change the limits of the sum over k. 
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The generating function for this restricted partition function can be obtained by omitting the factors (1 for 

j = 1,... in the product in the previous section for Q{n] X, 6)’s generating function or, equivalently, by dividing 
that product by these factors: 


N-l N-1 

X! 2/^ X! X,N -1 I qi,..., 

X—0 n=0 


N-l 

m—O 


1 

(1 + • • • (1 + yui‘ip ) 


Lemma 1 


Let uj = and let {gi, ■ ■ ■ ,qp \ 1 < qi < N — 1} he a set of p distinct integers. Then 


N-l 


P 


T.(-yr n 


1 


N—l—p 


m—0 


s=l m—0 0<Ki,...,Kp<m 


For the proof of this lemma, and that of Lemma 2 below, see Appendices C and D, respectively. 
Lemma 1 allows us to write eq.(lOa) as. 


Co.. 


•0 1- • -1 Ai ■ 




^ ^ - (ki + • • • + Kp)]. 


qi^qi 0<Ki,...,Kp<Mi 


(10b) 


We plan to use eq.(7) to perform the summations over the q variables in eq.(lOb). This is complicated by the fact 
that the sum in (7) is over the complete range of q, from 1 to TV — 1, but each g-sum above, while having this range, 
is restricted by the condition that qt qi for all i ^ 1. We will therefore express the restricted, p-dimensional sum in 
eq.(lOb) by a sum over unrestricted, lower-dimensional g-sums. An elementary example of this is 


1 


E Hx,y) = \ 


x<.y \ x,y X ) 

for a symmetric function h{x,y). Lemma 2 below generalizes this example to higher dimensions. 

We first define the p-dimensional vector : For each part Zk of a partition Z of p, with j parts, the p-dimensional 
vector Ufc has components 


1 if zi + ■ ■ ■ Zk i < zi + ■ ■ ■ + Zk+i i 


Uk = 


0 otherwise. 

Then for variables (gi, g 2 ,..., g^), is defined as 

Qz = liui H- qjUj = (gi,...,gi,g 2 ,...,g 2 , . ,qj,...,qj). 


Zl 


Z2 


q^ thus lies on a j-dimensional sublattice determined by the partition Z. 

Lemma 2 

Let h{qi ,..., qp) = h{q) be a completely symmetric function of the variables gi,..., g^. Then 

E ^(Qz)- 

l<qi<---<qp<M ^ ZGVp l<gi 
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The function f{qi, ■. ■, Qp) that appears in eq.(lOb) above, 

is however not a symmetric function on account of the factor in front. To correct this we divide the sum over qi,... ,qp 
in eq.(lOb) into pi subregions: 

f{qi,...,qp) = i X! + X! ^-^ f f{qi,...,qp), 

qi^qi [gi< 92 <'"< 9 p q 2 <qi<---<qp qp< -<qi J 

where the sums inside the brackets correspond to the p! orderings of (qi,... ,qp). Relabeling the summation indices, 
this is 

f{qi,---,qp) = X! {/(9i>■ • ■ H— + fiqp,-■■,qi) } = p'- fsiqi,---,qp) 

qi^qi qi<---<qp qi<---<qp 

where fs{qi, ■. ■, qp) is the symmetrized function, 

fsiqi,---,qp) = ^ 


(y^Sp 


Now, for (7 € Sp, define the p x p matrix 

Then 

(tCS'p 0<ki ,...,Kp<Mi 

where A = (Ai, A 2 ,..., Ap). Now defining as the integer in [0, A — 1] such that 
Ai'") = -(A • C>(cr)u,) mod N 

= ~ (^cr(ds + l) + • • ■ + ^cr((is+ 0 s)) Hiod A, with ds = Zl + ■ ■ ■ + Zg-l, 
we have, applying Lemma 2, 


.f{qu---,qp), = Y 


p+j 


{p;ki,k2, ■ ■ ■ ,kp)* ^ ^ 


p\ 


Y Y 


X Y Ml - (/ti + •••+ Kp) ]. 

0<Ki.. ,Kp<Mi 

The p sums over the k variables can be reduced to j sums over a variables, defined as 

ai = Kj H-+ 

Oik — _+ 1 + ■ ■ ■ + _’ 


aj = «p-.j +--- + S; 

Oi + • • • + CXj = Ki + • • • + Kp 


Then 


^(«l + -+«n)<?i+-+(«P-,-+- + «^p)92£f[ Ml - (ki + ■ • • + Kp) ] 

0</sii. .,«:p <Afi 

i / a.s + Z.S - I 

Y w“^«^+-+“^'«^A[Mi-(ai + --- + a,)]]^ 

0<ai,...,aj<Mi s—1 \ Zg 1 
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OCg -\- Zs 1 

where the multiplicative factor | | is the number of compositions of ag consisting of Zg or fewer parts. 

Zg-l 

Making this substitution, we can finally use eq.(7) to evaluate the sums over the g’s: 

(p; fci,fc2,.. .,kpy 


X! /(9i>--->9p), = (-If XI X 


(T^Sp Z^'Pp 

3 I ag + Zg-l 


p\ 

3 


^ X H[Mi-{ai + 

0<ai ,...,ocj<Mi 


n 


2c — 1 






We expand the product ni^i( ^ ~ powers of (—A^) as 


3 

r 

S=1 Ai,...,Aj=0.1 




(o-) ) — 


X 




and get, after rearranging the order of some of the sums, 

X = (-If X ^ ^ (_ivf.+-+A. 


Qi^Qj 


(T^S'P Z^'Pp 


pi 


3 / ag + Zg - 1 

^3) ] H 


■ 




X Y1 H|Mi-(o,+. 

0<ai,...,aj<Mi s=l V Zg — 1 / s=l 

Now if, for fixed (Ai,..., Xj), a particular At is equal to 1, then that term will contain a factor of d and we 

can replace the at in the binomial coefficient by . The above equation can therefore be written as 

(p; ki,k2,...,kp)* 


X /(9i>--->9p)> = (-1) 

Qi¥^Qj 


PHo 


X 


p! 


ZG'Pp 

X ^ ^ II[A4i — (ofi + • • • + Cij)\ 

0<Q:i. .,Q:j <Mi 


‘Pj_ 

Ai,...,Aj^0,l 


n 


a, + 2, - 1 \ ^ + 2s - 1 




2s - 1 


2s — 1 


Again take (Ai, ..., Xj) as fixed and assume that a particular At is equal to 1. The summand in the sum over the 
corresponding at is then H[M — (ai + • • • + C(j)]S^^ ^(<t) times factors independent of at- We can thus replace at in 

the step function by Xt'^K The sum over at is from 0 to Mi, and the sum over will be zero if > -^i- 

But, with the above replacement, the step function is also zero in this case. So the “zero condition” is enforced by H 
and we can simply make the replacement: H[M — {ai + ■ • • + aj)]5^ —)■ H[M — {ai + • • • + ' 

After performing the summations over all the cTs’s for which As = 1, we are left with the sums over the a’s for 
which As = 0. For Xt = 0, 


Ml / at + zt -1 

E 

at=0 \ Zt-1 


Ml — at — X! (1 — ^q)^q — X! ^q^q 

9=1 9=1 


M[ / at + Zt -1 

1 

at=0 y 2t — 1 

M( + (1 — At)2t \ 

{1 — Xt)zt j 

where the prime on the first sum over q means that the sum is over all q ^ t, and where 


M{ = Ml - X (1 - A9)«9 - X 

9=1 9=1 
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Summing again for As = 0, with M" = M[ — as, we have 


-^1 / cTs + ~ 1 \ / M” — as + (1 — Xt)zt 


=0 \ Zs-l 


H[M'{ - as] 


s - 1 / \ (1 - Xt)zt ) 

/ Qfs + 2s — 1 \ / — as + (1 — Xt)zt 

a,=0 \ 2s - 1 / \ (1 - Xt)Zt 


Ml + (1 — Xt)zt + (1 — As) 2 s 
(1 ~ Xt)zt + (1 — As) 2 s 




using the identity 


n m \ n — m 
m—0 \ I j \ h — I 


, 0 < I < k < n. 


Continuing, the effect of summing over these at’s is to remove a (—a) from the step function and the binomial 
coefficient, and to add (1 — Xt)zt to both the upper and lower elements of the binomial coefficient. So, after summing 
over all a variables, we are left with the factor 


H Mi-^XtX[ 


Ml + ELi(1 - ^t)zt - ELi \ ^ + 2s - 1 


ELi(l-AtK 


n 

s=l \ 2s - 1 


Putting this into the equation above for the sum over / and then that into (10b), and then making the replacement 
X](l - Xt)zt = p-E ^tzt, we arrive at 


_ (—1)^ (p; fci, fc2, ■ ■ ■, kp)* 

O-Ol-lAi-ApAp+i - ^ ^ p! 


^ (_jV)Ai+-H 


(y^Sn Z^T^p 


Ai,...,Aj=0,l 


r ^ j/A'-Mo-i-E's=iAA^i"^ + ^«)\^/xi"^ + 2s-r 
L ‘=1 J V Mi-E V -1 . 

For the terms for which Ai = • • • = Aj = 0, the sum over Z € Vp, and then trivially the sum over a G Sp, can be done 
using the identity 


^ ] (p! ki, ^2) ■ • ■; kp) — ^ ^ ] ^ki+2k2 -\—hpfcpiP ^ki+k2A —(Pi ki, k2, ■ ■ ■, kp) — p!; 

ZGVp j—0 fci,...,fcp 

(this is eq.(E4) in Appendix E). Then 


^0 - • - 0 1 • • • 1 Ai" • ApApj^-x 


{-1)N-Mo-^N (iV-Mo-1)! {p-ki,k2,...,kpy 

M2\---Mn-i\ 1 Ml! ^ A p! 

^ e ! p e o p 


Y^{-N)^ ^ H Mi-Y^XtXE^ 

/J-—1 Ai,..., Aj —0,1 _ t—1 

(N-Mo-l-Ei=iE{xi'^^+Zs)\ 3 /xi") + 2s-l\"^ 


Mi-E EX) 


n 

S=1 \ 2s - 1 


We now note that in the sum over Sp, the set of xi^^’s is invariant for all a that map a partition onto itself, 
ct( 0) = 0, since in this case a either just reorders the elements in each 0 € 0 or interchanges the order of the 0’s, 
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and neither of these operations change the set The sum over Sp and Vp can then be replaced by a 

sum over the set of partitions of the index set {12 • • -p} by multiplying each term by Zilz 2 '- • ■ • zjl ki\k 2 \ • ■ • kpl: 


E E - E ; 

zeVp aeSp ee-p{i2 -p} 

(p; ki,k2,...,kp)* ^ (p; ki,k2,... ,kp)* 

p! {p;zi,...,Zj) 




\ki 2^2 ■ ■ ■ pk 


With this substitution, 


Xi") ^ X{9s) = - 


^ I 

E I inod N, for ds = { , n^z} }■> 

, (o) 


and we get 


Co ---0 l - l,4i -,4pAp+i — 


(_l)^-^o-iAf I _Mo- 1)! 


M2\---Mn-i\ I Ml! 

Ai=l Ai,...,A3=0,1 


+ E (^1 - - 1)! 

eG-p{i2---p} 


Ml 




N -Mo-i-Y:>t^^\t{x{et) + zt)\ 1 ( x{es) + zs-i 


Ml - ELi ^tXiSt) 


n 


Zs - 1 


(lOd) 


Expression (lOd) above is an alternative expression for the coefficient to that stated in the theorem. The final 
expression is derived from (lOd) by the replacement of the sum over 7^(12 ■ ■ ■ p} by a sum over T^jAli • ■ • Ap'\, (which 
is a smaller set when the coefficients Ai,..., Ap are not all distinct), and then noting that if, for example, Ai = A 2 , 
a partition of the form {Ai ■ ■ ■ ){A 2 •••)•••(•••) will contribute the same to the sum as the partition obtained by 
interchanging 1 and 2. So, replacing the sum over 7^(12 • • -p} by a sum over the smaller set T^jAli • • • Ap} requires 
that the sum be multiplied by a factor of 2, or, more generally, by the factor M^J • • • M^^l. However, for partitions 
of the form {A 1 A 2 •••)•••(•••) or of the form {AiB ■ ■ ■ C){A 2 B ••• C) •••(•••), interchanging the indices 1 and 2 does 
not generate additional partitions if Ai = A 2 . To correct for this overcounting, the sum then needs to be divided by 
the factorials of the ng and the mo ^ multiplicities. A 


Corollary 5 // Oq + Oi + • • • + a^-i = 0 mod TV, then C[a] = 0 only if N divides ^ 


The proof follows by noting that every term inside the brackets in expression (lOd) is an integer, and that all terms 
other than the first one are proportional to TV. 


C. Coefficients of the form Co- -o i- 1 o- -ab 

In the special case where Ai = A 2 = ■■■ = Ap = a the expression for the coefficients simplifies considerably. In this 
case, in eq.(lOc), X^^ = {—aZg) mod TV, independent of the permutation cr. The sum over Sp can then be 

done and, after some manipulation, so can the sum over Vp. 
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Summing over the permutations we have, for Mi, = 0,1, 


Co---0 l - l a-'-ah — 


I (iV-Mo-1)! 


MJMbl 


Ml! 


zer„ 


M=1 Ai,...,Ai=0,l 


+ E iP',ki,k2,...,kp) 

3 

Ml - ^ AiX, 




N — Mq — 1 — X]t = l + Zt) \ 1 + Zs — 1 


A, ^ 


Ml - ELi 


n 


Zs - I 


We now replace the sum over (Ai, ... ,Xj) by a sum over (/3i,..., /Sp), which we define as follows: Let the part 
have multiplicity k (= kz„). Then, (after a possible relabeling of the parts), z„ = Zn+i = ■ ■ ■ = Zn+k-i- For fixed 
(Ai,..., Aj), let P equal the number of non-zero A’s in the subset (A„,..., A„+fc_i); i.e., /3 = E't^n ^ /3z„ is thus 

defined for each unique part z„. This definition can be extended to all integers q € [1,..., p] by setting Pq = 0 ii kg = 0. 


For q € [1,... ,p\, q and Xq will each appear Pq times in the sums over t, and the binomial coefficient 

will occur Pq times in the product over s. For a given (/3i,..., Pp), there are ^ ^ different (Ai,..., Aj)’s 

that correspond to it. Replacing the sums over (Ai,..., Xj) in the above expression by sums over (Pi ,..., Pp), 




P P 


^2^tZt, '^^XtXz-t — PtXt, 


t=i t=i 


t=i t=i 


and including this factor gives us 


C'o---0 l - la-.-ofa — (“1) 


N-Mo-l 


N I (A^-Mo-1)1 


MJMb! Ml! 


E E {P'^ki,k2,...,kp)* 

i=0 ki,...,kp 


ki 


/3i,...,/3p V Pi 


H 


Pp 


Mi-j^PtXt 


t=i 


N-Mo-l-ELif3t(Xt + t)\ P X, + s-l 


Mi-EEiPtXt 


n 


S — 1 


We now interchange the sums over the fc’s and the P's and perform the sum over the k's by using the lemma below: 


Lemma 3 For p> pi + 2p2 H--f pPp, 


. ki 

^ ^ y, 

j=0 ki,...,kp \ Pi 


(p;ki,k2,...,kp)* = 


pi 


1/5i/3i! 2/52/32! ••• p^r/3p!- 
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(For the proof, see Appendix E.) Then our result is 

N 


Cq-.-O l -la -ab 
Co-.-O l - l ha-'-a 


= (-1) 


N-Mo-1 


Ma + Mh 


Ma + Mh 

M„ 


N - Mo-1 

Ml 


n=l 0</3i,...,/3p<p 


H 


Mi-^PtXt 




N-Mo-l- ELi Pt{Xt + t) 

Ml — PtXt 

Xs = —sa mod N such that 0 < < A^ — 1; 

p = N -Mo-Ml-1. 


n 




Xg + 5 — 1 
5—1 


for = 0 or 1. 


(lOe) 


As indicated, this expression is valid for b < a as well as for b > a. The reason is that, in the proof of Theorem 3, 
in summing over the equivalence classes we have the freedom to “hide” any particular index in the ctq position, where 
it’s value has no effect, (other than on the requirement that it and the other indices satisfy condition (8a)). We can 
therefore sum over equivalence classes of the form [6, cti, ..., gn-i] for b either less than or greater than a. 


D. Coefficients of the form Co- oi- 1A1A2A3 

As previously noted, the two coefficients Cooiaas and CooDOiiiseS; corresponding to A = 6 and 10 respectively, 
are zero despite satisfying condition (8a). Generalizing, we consider coefficients of the form C'o...oi...iAiA 2A3 j with 
Mo, Ml > 1 and 2 < Ai < A2 < A3 < A — 1; i.e., coefficients with non-zero Mq and Mi that contain 3 and only 
3 indices greater than 1. We are interested in the question of when coefficients of this form equal zero even if they 
satisfy conditions (8). From (8b,c) we have 

Ml jAi -f A2 “t“ A3 = 0 mod X, 

Mo + Ml+3 = N, 
and from eq.(lOd), 

C0...01...1A1A2A3 tx { (Mi + 2)(Mi + l)-Ai 7 [Mi-A(AiA 2 )](A(AiA 2 ) + l) 

-N{H[Mi- X{Ai) ] ( Ml + 1 - A(Ai) ) + (Ai ^ A2) ) + N^H[ Mi - X{Ai) - ^(Aa) ] 

For Co-.-oi---! A1A2A3 to be zero, N must divide (Mi -|- 2)(Mi -|- 1) for Mi < N — 4. If A divides (Mi -|- 2)(Mi -|- 1), 
then it also divides (Mq -|- 2) {Mo + 1): 

(Mi+2)(Mi + l) (iV_Mo-l)(A-Mo-2) ^ , (Mo + 2)(Mo + l) 

- N -=- N -= A-2Mo-3+- - - . 

The X functions in eq.(ll) are: 

A(Ai_ 2 ) = A — Ai^ 2, 

r Mi-fA 3 ifA 3 <A-Mi; 

A(AiA 2) = —(Ai -I- A2) = Ml -I- A3 = < 

[ M 1 +A 3 -N if A3 > A - Ml. 

If A3 (and therefore Ai and A2) is less than N — Mi, all of the Heaviside step functions in eq.(II) are zero and the 
coefficient is therefore non-zero. So we will assume that A 3 > N — Mi. That leaves only 4 possibilities: 

(1) N — Ml ^ A3 and Ai ^ A2 ^ N — Mi', 

(2) X — Ml ^ A2 ^ A3, Ai < X — Ml, and Ai -t- A2 ^ 2X — Mi', 

(3) X — Ml ^ Ai ^ A2 ^ A3 and Ai -t- A2 ^ 2X — Mi', 

(4) X — Ml ^ Ai ^ A2 ^ A3 and Ai -t- A2 ^ 2X — Mi. 
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Now let ^ be the symmetry operator defined as multiplication of an index set [a] by (TV — 1) followed by addition 
of 1. Under this operator, 

Mq —^ Mq = Ml, Ml — y Ml = Mq, Ai — y A'l = N — A^ + 1, A2 — y A2 = N — A2 + 1, ^3 — y A3 = N — Ai + 1. 

Thus, if TV — Ml < Ai < A 2 < A 3 , then A^ < A 2 < Ag < From the discussion above and from Proposition 

2 , in this case, 

Co. A2A3 (-i)'^-'c^ *[0---01---1 AiA 2A3] 0i 

and we need to consider only cases (1) and (2). Setting the right side of eq.(ll) to zero in these two cases gives us 
the corollary to Theorem 3: 


Corollary 6 For N = Mq + Mi + 3 with Mq, Mi > 1, the coefficient Cq... 01...1 A1A2A3 satisfies conditions (8) but is 
zero if N divides {Mi + 2)(Mi + 1) and either 


Ai < A 2 < N - Ml, Ai + A 2 = TV + 1 - (^1 + 2)(^i + 1) ^ ^3 = Mo + 2 + 


N 


TV 


N - Ml < A 2 < Aq, A 2 + A 3 = tv + 1 + (^0 + 2)(^o + 1) ^ Ai=Mo + 2-^°^ 


N 


TV 


By Proposition 2, all coefficients that are related to these zero coefficients by an additive- and/or a multiplicative- 
symmetry operator also equal zero. 

For integers a and q, TV in the corollary has the form TV = q{q -|- l)/a, with values 6,10,12,14,15,18,..., corre¬ 
sponding to sequence A139799 in OEIS Q- 

As an example, we take TV = 6. In this case TV divides (Mi -|- 2)(Mi -|- 1) = 6,12, 30,42,..., but only the first two 
integers satisfy the condition TV > Mi -1- 4. From the corollary, 

C001245, C011235 (= —C^*[ooi 245 ]); 

C 001335 , C 011244 (= —C^*[001335]) 

are zero coefficients. Including all the coefficients generated from these by the Proposition-2 symmetries, there are a 
total of 12 zero coefficients. 


As a second example, TV = 10. TV divides (Mi -|- 2) (Mi -|-1) 


Coo00111377j Cooo1111448 
Coo00111368) Cooollll358 
Coo00111458: Cooollll367 


(= ~C^*[ooooill377])) 
(= ~C^*[ooooill368])) 
(= —C^*[0000111458])- 


20, 30, and the corollary gives us the zero coefficients 


Including the Proposition-2-operator-generated coefficients, there are 120 zero coefficients of this form. 


IV. ADDITIONAL RESULTS: MULTIPLET STRUCTURE 


Having a formula to calculate any coefficient in the expansion (3a/b), it is of interest to determine how many 
calculations one needs to perform to solve for the complete determinant. The number of solutions to eq.( 8 b,c) is given 
by the expression [^, @ , 




d\N 


N 



(A003239)Ar+i in OEIS. 


( 12 ) 


where (j) is the Euler totient function Q- Then, except for coefficients that equal zero as a result of satisfying the 
criteria of Corollary 6 and any other “accidental zeros”, eq. (12) gives the number of terms in the expansions (3a/b). 
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As discussed above following Proposition 2, coefficients can be grouped into additive multiplets, related by the 
additive symmetry, and into super-multiplets, related by the additive symmetry and/or the multiplicative symmetry. 
For example, for N = 5, the coefficients can be formed into the table: 


5 : 

{C'sooools 

311 : 

{^siooils) {C'aoiiols 

221 : 

{C'i2002}5) {1^10220 Is 

mil : 

{Qllllll- 


The set of coefficients consist of six additive multiplets (where n is the number of elements), and four 

super-multiplets. (In this example there is only one super-multiplet for each of the allowed partitions of 5, but that 
is not generally the case; see the table for the N = 8 coefhcients in Appendix G.) 

The nonnegative integers less than N form a group under addition modulo N, which we will denote as (Z/A^Z)+. 
The group (Z/A^Z)+ is of order N, and is a cyclic group, so all its subgroups are cyclic. The set of integers less than 
and relatively prime to N is the multiplicative group of integers modulo N, denoted as (Z/NZ)^. This group has 
order (j){N). The direct product of these two groups, 

Hn = iZ/NZ)+ X (Z/NZ)^, 
is thus of order N(I){N). 

If the index set [M] is such that a * [M] = [M] for all elements a in a group, we say that [M] is invariant 
under the action of that group. To find the numbers of additive multiplets and of super-multiplets, we will need 
to count the number of index sets that are invariant under the various subgroups of (Z/A^Z)+ and of H^, respectively. 


A. The number of additive multiplets 

From eq.(3b) the determinants are expressible in the form 

det [xo,a;i,...,a;iv-i] = ^ • 

For a given coefficient , the elements of the n-element additive multiplet are obtained 

by circularly permuting the index set [MqMi ... 




Mq M\ • • • Afjv — 


= { c 


Mq Ml • • • Afjv 


I 

^P[MoMi---Mn-i]^ ^P'^[MoMi---Mn-i] ' - ' ^^P'^-'^[MoMi---Mn-i] } 


where P is the circular permutation operator defined in eq.(9). We define the symbol [ ] as equal 

to the value of the coefficient that labels the additive multiplet: 


[ {C'moMi-'-Mn.i}" ] — 

Then, as a sum over these multiplets, the above expansion can be written as 


(13) 


det [a;o,Xi,...,a;7v-i] = ^ [ {C'MoMi.-Mw.Jn ] 

Multiplets 






V + (-i)^ 




Mn-1 

^0 


+ ••• + (- 1 ) 




Mn_ 


(14) 


Proposition 3 For a given N, the number giq{n) of n-element additive multiplets is 

N-n / 2fi 


9Nin) = 


l + (-l) 


4n^ 


i:(-i)-v(5) 


d\n 
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Proof: 

Operating on an index set [M], the additive group (Z/7VZ)+ generates a set with N not-necessarily-distinct elements: 

{Z/NZ)+ *[M] = { I 0 < a < iV - 1}. 

Let {M}n = { [M], P[M], P‘^[M],..., } be an n-element additive multiplet. Then n divides N and, for 

d = N/n, each element in the multiplet is invariant under the action of any element of the subgroup S'j’ C Z/NZ^. 
One of the index sets in {M}, which, without loss of generality we can take to be the “labeling” set [M], must then 
have the form 

[Mq, Ml ,..., Mjv_i] = [Mo, Ml,..., M„_i, Mq, Mi, ..., M„_i, • • • , Mq, Mi, • • • , M„_i], 
in which the subset {Mq, Mi, ..., M„_i} is sequentially repeated d times. We have then: 

N-l n-1 

^ M, = {N/n) ^ M„ 


qMq = Y q+iq + n) + iq + 2n)-\ - \-{q + n{ N/n - 1) ) ] 

^71 9 


= - y. q^q ”1— 

Tt ( ^ 71 


N N-n 


n 2 


YMq. 


Conditions (8b,c) then require that 

n —1 

YMq = n, (15a) 

qr=0 

and 

Y^ q^q + ^ ^ 

The second of these equations can only be satisfied if TV — n is even. (If fV — n is odd, n must then be even in order 
that the right side to be an integer. But then N is odd, which contradicts the requirement that n divide N.) In this 
case then, 

n—1 

^ gMg = 0 mod n, fV — n even. (15c) 


From eq.(12), the number of solutions [Mq, Mi, ..., M^-i] that satisfy both eqs.(8a,b) and (15a, c) is therefore 


FN{n) = 


l + (-l)^ 


E l f 


and n-member multiplets {M}„ and {C*j^-^}n exist only if n is a same-parity divisor of N. 

The number of solutions Fn {n) is a sum over all coefficients that lie in d-element multiplets for all d that divides n 
and is of the same parity as n and N. Let l{d) denote the total number of coefficients C[m] for which d is the smallest 
integer such that P‘^[M] = [M]. Then the number of n-element additive multiplets, (/Ar(n), is l{n)/n, and 
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where the notation d\\n means that the sum is over all same-parity divisors of n. It is straightforward to show that 

d\ |n 

by considering separately the n odd and the n even cases and using Mdbius inversion. Then 

n /n\ , f 


9N{n) = 


4n^ 


■1:1: 5'*© 

d\\n q\d 


It remains to prove the identity 

'd 

.9 


d\\n q\d h\n 


(16) 


where f{q) is an arbitrary function of q. In evaluating the right-hand side of this equation, we express n as n = 2^m, 
where m is an odd integer and k > 0. The divisors of n are of the form b = 2©, where a divides m and is therefore 
odd, and I <k. Then 


r /(a), k = Q] 

RHS = ^(-I)"+V f{h) = ^ M (^) J /(2a) + /(a), k = I; 

b\n a\m I 

I /(2©) -/(2'=-ia), k>2. 


(17) 


We now consider the expression on the left side of eq.(16), and first take n to be odd, {k=0). Then all of its divisors 
are also odd, and we can make the replacement fi||n —>■ d\n. For a given divisor 5, we set d = be, c < n/b, and isolate 
the q = b term in the sum over q: 


LHS, 


odd — 


d\n q\d 


El: /(")—+/<'» E 


= •••-k 


c\{n/b) 

(?) 


H-= 


/©) + •••. 


a—b\m—n 


where we’ve used the identity 

X! ~ ^ {~) = k-im), 

c\m 

which follows from Mobius inversion of 

Hm) (18) 

d|m 

Eq.(16) is therefore valid for n odd. 


For even n, we again set n = 2^m, with k > 1 and m odd. Its same-parity divisors are d = 2*a, with ajm and 
1 < I < k. Setting q = 2®p, with 0 < s < ? and p an odd integer that divides a, the left-hand side of (16) is 


2'’’-'m /2'=-' 


m 


2^-^a 


P 


/( 2 ©). 


LHSn even “ X/ X/ X/ X/ ^ 
a\m 1 — 1 p\a s—0 

But since p{2^m) = 0 for any K > 1, the only non-zero terms in the sum over d\\n are those for which I = k, fc— 1 > 0; 


LHS„ 


= EE 

a\m p\a 


E TTI (Vfl 

— ll\ — 
a \ a 


.s=0 


k-1 


/(2*p) + (l-di©^ — — 

a \ a 

s=0 ^ 


1—s, 


/( 2 ©) 
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The two identities 

fx{2x) = —^{x), 4>{2^x) = 2^~^(j){x) for K > 0, 

where x is in both cases an odd positive integer, can be used to take the Mbbius and Euler functions outside the 
parentheses: 




a\m 


p\a 


k-1 


k-2 


/(2V) + E - (1 - Kk) 2/(2'=-V) + E 2'=-^-V(2^p) 


s=0 


s=0 


The two sums from s = 0 to A: — 2 cancel and we have 
LHS n eve 


aim p|a 

Now applying the “odd-n” result, 

d\n q\d 


/ [ /(2"p) + (2<5i.fc - l)/(2'=-ip) ] . 


EE ^(5)^9 = 

d\n q\d b\n 


to the last expression, it becomes 


LHS 


n even 


(^) [/(2'=a) + (2<5i,fc - l)/(2'=-ia) ], 

a|m 


in agreement with the 2nd and 3rd lines in eq.(17). A 


Note that the expression 



corresponds to the sequence A131868 in OEIS. So gN{n) = (A131868)„||Ar, and the total number of additive multiplets, 
riAM, is 

J29N{n) = ^(A131868)„. 

n||Af n||Af 


B. The number of super-multiplets 

To find the number of super-multiplets, we consider the direct product: 

Hn = {Z/NZ)+ X (Z/NZ)^. 

We will denote an element oi either with boldface, a, or as the pair (a, b). Group multiplication of two elements 
ot Hn is defined by 

{a,b) ■ {c,d) = {a + c,bd); a,c€(Z/AZ)+, 6,de(Z/AZ)^. 

In this formalism, the symmetry operations in Proposition 2 are expressed as the action of an element (a, b) on an 
index set [M], which we will denote as (a, b) * [M]: 

(a, b) * [M] = [Mjv-a, M 2 b-i-a, ■ ■ ■ , Af(jv-l)&-i-a]- 

In particular, if & = 1, 

(a,l)*[M]=F-“[M]. 
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A direct-product group Gi x Gi is cyclic iff Gi and G 2 are both cyclic and have orders that are relatively prime 
to each other. By Gauss’ Theorem, (Ti/N'L)^ is a cyclic group iff N is of the form 2, 4, p" or 2p" where p is an odd 
prime and n > 1 Q. Consequently Hfq is cyclic only for the case where is a prime number. 

We will denote an order-d subgroup of (Z/AfZ)+ as iSj"; since (Z/A^Z)+ is a cyclic group, all its subgroups are 
cyclic, and each d corresponds to only one subgroup. Likewise, we will denote order-d subgroups of (Z/A^Z)^ as S^.^, 
and those of idjv as Sd-^a, where the a indices are used to distinguish same-d subgroups from one another. 

Operating on an index set [M], idjv generates a set with N(f){N) elements: 

* [M] = { (a, 6) * [M] I 0 < a < AT - 1; 1 < 6 < Af - 1, gcd(6, A^) = 1 }. 

If [M] is a member of an m-element super-multiplet, then [M] is invariant under the action of an OTdeT-(N(f>{N)/m) 
subgroup of H^. If [M] is invariant under the action of all elements of a subgroup Sd-a, then cj- * [M] = Cj * [M] 
for all elements c^, Cf, of a coset {cfc} of Sdia- If the elements of the set { [M], Ci * [M], ...,Cm-i * [M] } are 
all distinct, where m = N(j){N)/d and cj- is any element of the coset {c^}, then this set is an m-element super-multiplet. 

For Sd-a Q Hn we define K[Sd-,a) to be the number of index sets [M] that are invariant under Sd-,a and in addition 
satisfy eqs.(8b,c): 


K{Sd-,a) = #{ [M]\ Sd-,a * [M] = [M], C Hn] Y,kMk = 0 mod N- Y,Mk=N 

{ k k 

If [M] is invariant under the action of the subgroup Sd-,a, it is also invariant under all subgroups of Sd-^a- It may also 
be invariant under larger subgroups of Hn that contain Sd-,a', i-e., it may have a higher degree of symmetry than that 
required to be invariant under just Sd-,a- K{Sd-,a) therefore counts these higher-symmetry sets as well. We define 
L{Sq-^p) as the number of index sets satisfying (8b,c) for which Sq-^p is the largest subgroup that they are invariant 
under. Then 

K{Sd-,a) = Y. 


Let 


N(j){N) 
n = - 

q 


Then L(Sq)/n is the number of n-element super-multiplets, and the total number, nsMy is 


nsM 


E 

SqCHN 


qHSg) 


(19) 


We can formally solve this equation by expressing the previous equation in matrix form: 

K = AL, L = A-^K, 

where the elements of the matrix A are either 1 or 0 depending on whether or not a particular subgroup Sq-^p contains 
the subgroup Sd-,a or not, and the rest of the notation should be clear from the context. 

Because of the more complicated structure of the multiplicative group, calculating the number of super-multiplets 
is more involved than for the additive multiplets. Expression (19) however simplifies in the cases where N = p ov 2p, 
where p is an odd prime number. In Appendix A, we calculate the number of super-multiplets for these two cases. 
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Our results are: 


nsM{N =p) = 

nsM{N = 2p) = 


1 J 

f 

1/ 

^2p\ 

p{p “ 1) 1 

2p 1 

\ P ) 

1 

J 


/ 4p \ 

2p{p - 1) 

^ 4p 

\2p) 


p^ — 1 

- -^ 

P 

Ap^ + 1 
Ap 


m<p—1 

^ E 


p-i 


2p \ o /'„2 


2m 


m 


2(p2 _ 1) 


p 


p 


m<(p—l)j2 

+p E ' 

m\{p-l) 


p-1 


m 


p + Am +1 1 — (—!)■ 


2m + 1 


Using eq.(12) and another series expression in OEIS, 
2m 


1 A ( ^ \ 

2 \m/ 


= (A081875)„, 


m\n 

expression (20a) becomes 


4m 

2m 


1 


2p-2 


nsM{N = p) = , " ' { (A003239)p+i + 2p(A081875)p_i + p-l-p 

PyP- [ \ P-I 


(20a) 

(20b) 


l-(-l) 


ezlL ( 2m 


m 


V. CONCLUSION 

The expression for Cao---aN-i in Theorem 3 constitutes in principle a complete solution to the problem of calculating 
the determinant of an arbitrary circulant matrix. This expression has some similarities to that of Wyn-jones, as well 
as some significant differences, but the proofs of the two expressions are along completely different lines. 

However, Wyn-jones’s method does require that one find all of the null subsets/multsets contained in a generally 
much larger set than the set of partitions employed in Theorem 3. And since, in applying Theorem 3, one deals only 
with partitions of sets that do not include any a = 1 indices, it may be that it is more efficient than the method of 
Wyn-jones for coefficients with Mi >> 1. 

In many cases, the number of terms in the expression for a particular coefficient can be decreased by using one or 
both of the symmetry operations in Proposition 2 to equate the coefficient (up to a sign) to another coefficient. This 
can often be done in one of two ways: By either making the new multiplicities Mq and/or as large as possible, 
thereby making the set of partitions to be summed over smaller, or by setting to 0, making the constraint 
Ml — ^ \sX{9s) > 0 more restrictive. 

To illustrate this second method, let C'oo...on_i be such that the integer n is not in the index set. Then if we add 
(iV -I- 1 — n) to [oo,..., oat-i], the resulting index set will not contain 1, so M[ = 0. We can always do this except for 
the index set [0,1, 2, • • • , iV — 1]. 

Consider the previous example, CooiiiisTss- The missing indices are n = 2,4, 5, 6,9, so TV -I- 1 — n = 9, 7,6, 5, 2, and 

Cooilll3788 = —^0000267799 = —<^0455778888 = +<^3446677779 = —<^2335566668 = +<^0022333359- 

Of these, <70000267799 is the easiest to compute as there are only two partitions, (677)(29) and (677)(2)(9), that 
contribute to the sum. 

This method of setting Mi to zero has some similarity to Wyn-jones’ method, in its reliance on the null sub¬ 
sets/multisets. The contributing partitions of the set {267799} would be in that case the trivial partition, (267799), 
and the partition (677) (299). 
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For CooiiiisTsSj the number of non-zero terms in the sum can actually be reduced further by instead multiplying 
the index set by {N — 1) and then adding 1. This has the effect of interchanging Mq and M^: 

[0011113788]''® ^ [0022379999]+^ ^ [0000113348] 

In this case the only contribution to the sum comes from the trivial partition (334), with X(334) = 0. All other 
partitions correspond to values of X{9) that are greater than M[ = 2 and therefore do not contribute. 


Appendix A: Derivation of eqnations (20) 


We will use the notation 

(5l/ = l)^{l,5,5^•••,ff"-' }, 

or the short-hand notation {g\g‘^), to denote a general order-d cyclic group generated by g. 
Throughout this section, p will denote a prime number greater than 2. 


Proposition 4 For N = p or N = 2p , 

The proof of this proposition uses the following lemma: 

Lemma 4 .• Let f and g be mappings from the set of subgroups {5'c;} of a cyclic group G into the set of nonnegative 
integers Z such that 

fiS,) = Y. 9iS,)- 


Then 



In both equations, the sums on the right are over all subgroups that contain the subgroup on the left. 

Proof of Lemma 4: Let G be of order k. Since G is cyclic, for each integer d that divides k there is one 
and only one subgroup of order d in G. Then we can make the replacements f{Sd) —>■ fid), giSq) —>■ giq), where 
fid), gid) are defined for all integers that divide k. Then the condition in the lemma becomes 

fid)\d\k = Y 5 ( 9 ) 

d\q\k 

where the sum is over all q that divides k and is divisible by d. For a given d and q that divide k let 



and 

fim) = f , gin) = g . 
Then the sum above becomes 

Jim) = E gin), for m\k. 

n\m 
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Now let /(s) and g{s) be extensions of f{s) and g{s) to the full domain of the nonnegative integers which preserve 
this relation but which are otherwise arbitrary: 

/(s) = f{s) ) 

> for s|fc; 

?(s) = 9 is) J 

f{s) = '^g{a) for all s G Z. 

a|s 

By Mobius inversion, 

sis) = ^/r(s/a)/(a). 

a|5 

Setting s = n = k/q and a = m = k/d we get 

did) = t^Wd)fid) 

q\d 

and 

diSq) = 9id/q)fiSd)- 

Sd^Sg 

A 


Proof of Proposition 4 Let D = N<p{N). 

If iV is a prime number then Hn is cyclic and from Lemma 4 we have 



d\D q\d 


= i J2d^i^<i)^id)^ 

d\D 

where the last line follows from the eq. (18) identity. Since each divider d corresponds to one and only one subgroup 
Sd, this sum is equal to the sum in the proposition. 


Now let N be equal to twice an odd prime number. In this case TLat is not cyclic. We consider the element 
(1, b) € Hn, where b is any element of (Z/NZ)^, and let [M] be an index set that is invariant under (1, b). We have 


[Mo, Ml, M 2 , ..., Mjv-i] = (1, b) * [M] 


[M^r-i, Mt,-i_i, M2f,-i_i,..., M(^_i)f,-i_i], 


and so Mq = Mn-i and Mk = for all A: = 1,..., iV — 1. & is an odd integer, so if k is even, {kb~^ — 1) is odd, 

and vice versa. Each index with an even subscript is thus equal to an index with an odd subscript, and we can write 
the sums in {8b, c) as sums over just the even indices: 


N-l 

N=YgMk 

N-l 

kMk 

k=0 


(N-2)/2 

2 y]] M 25 , 

q=0 

(N-2)/2 (N-2)/2 

Y, i2q + 2qb-^ -l)M2q = 2{l + bN Y dM2q 

q=0 q=0 


N 

Y’ 
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Since by assumption N is twice an odd prime, ^ kMk is an odd integer and is not congruent to 0 mod N. Thus, 
index sets that are invariant under any subgroup Sd that contains an element (1,6) for any b € (Z/A^Z)^, (which 
of course includes the full group H^), do not satisfy condition (86) and, as a a consequence, K{Sd) = 0 for such a 
subgroup. Then 

^ ct>{d)K{Sd)= ^ ct>{d)K{Sd) 

SdQHN SdQHN 

where the proper subgroup 

Hn = {0,2,A,...,N-2}x {1, 3,5,...N/2-2,N/2 + 2,...,N -1} C Hn 
is cyclic. The proof then follows as above for the N equal to a prime number case. A. 


If [M] is invariant under a cyclic subgroup Sd of Hn, then the conditions a* [M] = [M] for all a € Sd can be 
replaced by the single equation g* [M] = [M], where g is a generator of Sd- The product subgroup 5”^ x S'f is cyclic, 
with generator (1,1), since it can be identified with the cyclic group 5')() = (Z/iVZ)+. Its generator, expressed as a 
permutation on [M], is 

(1,1)* = (0123---(A- 1)) 

The only solution to (1,1) * [M] = [M] and to N = (eq.(8c)), is [M] = [111---1]. This solution however 

satisfies conditions (8b) only for odd N. Consequently, 

( 1 for odd A, 

K{S+xS^) = \ (Al) 

I 0 for even N. 

In addition, since all index sets are invariant under the identity subgroup Ai = S^ x S'f = {(0,1)}, 

We will use the lemma below to calculate nsniN) for N = p and N = 2p: 

Lemma 5 If [M] is invariant under the action of the group element {a,b), bf^\, then P'^[M] is invariant under the 
action of (a + (1 — 6“^)c, 6) 


Proof: Acting on the multiplicity index set [M] = [Mq, Mi, ..., the circular permutation operator permutes 

it to 


P"[M] = [Mn-c, Mi_„ ..., Mn-1-c] 

Then 


(a + (I- 6 -^ 0 , 6 )*(P‘=[M]) 


(a + (I — 6 ^)c, 6 ) * [Mn-c, Mi_c, ..., Mn-i-c] 

[]\d^N—c)b~^—a—c-\-b~^c: ^{l — c)b~^—a — c+b~^c7 ■ • ■ i -^(N—l — c)b~^—a—c+b~^c] 
[Af—a—C5 ^b~'^—a—ci • • • i ^(N—l)b~'^—a—c\ 

P‘=(a,6)*[M] =P"[M] A 


a. N =p 


For N = p, we have, from Proposition 4, 

nsMijp) = , , X! <P{k)K{Sk), 

Hp = {0,l,2,...,p-I}x{l,2,...,p-I}. 
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Hp is cyclic. All of its subgroups are therefore cyclic and are equal to the direct product of a subgroup of (Z/pZ)+ 
and of (Z/pZ)^. The subgroups of (Z/pZ)+ are = {0} and = {0,1, 2,... — 1}, with generators 0 and 1, 

respectively. The subgroups of Hp are then of the form 

S+ xS^ ={0}xS^ = {{0,g)\i0,gr), 

and 

={0,1,2,...,p-l}x =((l,g)|(l,gr), 

where d divides {p — 1) and g is a generator of C (Z/pZ)^. Then 


1 

nsMip) = -7-^ KiSt xS^) + ip- 1)K{S+ x S^) + ^ miKiSt xS^) + {p- l)KiS+ x S^)] 

(A3) 

where we’ve used the relation (j){pd) = (j){p)4>{d) = {p — l)(/)(d). The d=l functions, K{Si x S^) and K{Sp x S^), 
are given by eqs.(A2) and (Al), respectively, and we have 



K{St xS^) + (p- l)K{S; X S^) 



For d > 1, we show that 


K{S; xS^)= K{St X S^) 


2(P-1) 

d 


P-1 

d 


(A4) 


To prove the hrst identity, we set (a, &) = (0, g) and (a+ (1 — 6 ^)c, 6) = (1, p) in Lemma 5 above. Then c = 5 ((/ —1) 
For each [M] that satisfies {0,g) * [M] = [M], there is a circular permutation [M'] = [M] that satisfies 

(l,g) * [M'J = [M'J, and vice versa. The number of solutions to both eqs.(8b,c) and to (l,g) * [M] = [M\ is therefore 
the same as the number of solutions to (8b,c) and {0,g) * [M] = [M], and so K{Sp x S'^) = K{Si x S^). 


[Note that this argument fails for d = 1 since in this case 5 = 1.] 

To prove the second identity, we identify x with and set m = (p — l)/d. has (m — 1) cosets that 
we denote as {ci},..., {cm-i}, with Ck € {ck} C [l,2,...,(p— 1)]. Multiplication mod p by p induces a cyclic 
permutation g* on [012 • • • (p — 1)]: 

5 * = (0)( 15 • • • g ^-^)(Cl (pci) • • • ig‘^~^ci )) • • • (Cm-i (gcm-i) ■ ■ ■ )■ 

Then [M] is invariant under g* if and only if 

Ml = Mg = ■■■ = Mgd-1, 

= Mgei = • • • = (A5) 


Mcm-l — MgCm-1 — ■ ■ • — Mgd-li,^_^. 


Therefore, 

N-l 


^ ' Mq — Mq + d{Mi + Mci + • • • + 

q=0 

and 

N-l 

qMq = ( Ml + ciMci + • • • + )(1 + g + ■ • ■ + g^ ^). 

9=0 


g + ---+g 


d-l 


g^-1 

5-1 


kp 

5-1 


We have 
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for some integer k, where the last equality follows from g‘^ = 1 mod p. {g — 1) then divides kp; since p is prime, it 

must divide k. Therefore, for(i> 1, l + f/H-h g‘^~^ = 0 mod p and condition (8b) is satisfied for any index set [M] 

that satisfies equations (A5). 


Condition (8c) and the first equation above require that 


Ml + Me 


-Mr 


N- Mo Mo - 1 

= m — 


d d ’ 

so we must have that 

Mq = 1 “t“ deX ^ 1, Ml Mq^ “t“ * ' * “t“ ^Cm — l ~ ^ ^ — 0; 
for some integer a G [0, m]; (Mq = 0 is excluded since d > 1). For a given value of a, the set {Mi + 1, 


Mr 


1 } 


is a composition of {2m — a). The number of compositions of an integer n into k parts is ^ ^ | ^, so the number 

of solutions to the set of equations in (A5) which also satisfy (8c) is 


K{st y<sS) = f2 


a=0 



(A6) 


Putting equations (A4) and (A6) into (A3), we get eq.(20a). 


b. N =2p 


In this case, 

nsM{2p) = 2p(p- 1 ) 4>{k)K{Sk), 

H 2 p = {0,2,4,...,2p-2}x{l,3,5,...p-2,P + 2,...,2p-l}. 

H 2 p is cyclic and its subgroups are of the form 

St xS^ ={0}xS^ ={i0,g)\{0,gr), 

St xS^ ={0,2,4,...,2p-2}xS^ ={{2,g)\{2,g)P'^), 

where ={!,(/,... ,g‘^~^} C (Z/(2p)Z)^ is a d-element cyclic subgroup generated by the odd integer g. As indicated 
above, St x and 5”+ x are both cyclic, with generators (0,p) and (2,p), respectively. Then 


nsM{2p) = 


d>l 


2p{p - 1) 


K{St X St)+ {p- l)K{St X St) + Y. ^(d) [K{St X St)+ {p- l)K{St x | . 

(A7) 


d\{p-l) 


For d = 1 we have 


1 


A'(S.+ xSP = j;: 

^ (?=l,2,p,2p 


2q 

q 


2(p-l) , 1 

p 4p 


2p 

P 


4p 

2p 


The cyclic subgroup 5+ x St has the generator (2,1). As a permutation on [M], 

(2,1)* = (024 • • • 2p - 2)(135 • • • 2p - 1). 

The only solutions to (2,1) * [M] = [M] which satisfy (86, c) are [M] = [2020 •••20] and [0202 •••02], and so 
K{S+ X St ) = 2. 
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The identity K{Sp x S^) = K{Si x S^) ioi d > 1 follows by a similar argument as that in the previous section, 
using Lemma 5, with c = 2g{g — 1)~^. It remains then to calculate K{S^ x S^) for d > 1, which we do in essentially 
the same fashion as was done for the N = p case. We identify x with and consider the permutation 
( 0 , 5 )* = g* acting on [012 ... — Ij: 

9* = (0)(l5 • ■ ■ 5 • Cl • ■ • g'^~^ • ci) • ■ • (c^-i g ■ Cm-i ■ ■ ■ 9'^~^ ■ Cm-i) 

{2 2g--- 2g‘^~^){2ci 2gci ■ • • 2g‘^~^ci) ■ ■ • ( 2 cm-i 2gCm-i ■ ■ • 2g‘^~^Cm-i){p)- 

Again the Cfc’s label the cosets of and m is defined as {p — l)/d. If [M] is invariant under then g * [M] = [M], 

and conditions ( 8 b,c) become 


2p — 'y [ Mq — Mq + Mp + d[Mi + Mci + ■ ■ ■ + + M 2 + M 2 C 1 + ■ ■ ■ + M2cm-i) 


or 


Ml + ALci + • • • + + M 2 + M 2 C 1 + • • • + M2cm-i — 

and 


2p — Mq — Mp 


= 2m — 


Mo + Mp-2 


(A 8 a) 

(A 8 b) 

(A 8 c) 


d d 

0 = qMq = pMp + (1 + g + ■■■+ g^^ ^)( Mi + CiMcj + • • • + Cm-iMc^_^ 

+ 2 M 2 + 2ciM2ci + • • • + 2Cm-lM2cm-i )■ 

Since g‘^ = 1 mod 2p, g^ = q{2p) + 1 for some integer q. Then, for d > 1, 

T.(SP = l+s + ^- + /-' = ^ = ^- 

The trace of is a sum of d odd integers, which is odd or even depending on where d is odd or even. Since p is odd, 
the parity of the integer 2q/{g — I) is the same as the parity of d. As a result. 


Tr{S^) ^ 


0 mod 2p if d is even; 
p mod 2p if d > 1 is odd. 

We consider the cases d equal to an even integer and d equal to an odd integer greater than 1 separately: 


(I) d = 2k. Eq. (A 8 c) then requires that Mp is even. From (A 8 a/b) Mq must also be even. Setting Mq = 
2mo, Mp = 2mp, we have from (A 8 b) 

Ml + Mci + • • • + Mc^_i + M2 + M2C1 + ■ ■ ■ + M2cm-i ~ 2 m — a 

where a is an integer less than or equal to 2m defined as 

Mq + Mp — 2 mo + mp — I 
d ^ k ’ 

or mo + mp = ka + 1, > 0 . 

The set {Mi + I, Mc^ + l, ■ ■ ■, Af 2 c„,_i + 1} is then a composition of (4 to —a). For each a there are (fcQ; + 2) combinations 
of Mo and Mp such that Mo/2 + Mp/2 = ka + 1. For k > 1, a can take on values from 0 to 2m, while for /c = 1 the 
range of a extends downwards to -I. The total number of solutions {Mo, Mi,..., M 2 c^_j^} is then 


K{S+ X S^)\d even 


2m / 4 to — a — I \ 

'y^{ka + 2) j j + 5d,2 

a=o y 2m — 1 J 




(P- 1) 


4m + 1 
2m + 1 


4m 

2m — I 


d = 2 . 


d > 2 . 
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(II) d = 2k + 1, fc > 0. In this case, 

qMq = p{Mp + Ml + ciMci H-h ) mod 2p, 

and condition (8b) requires that the sum Mp + Mi + CiM^j H-h Cm-iMc^_^ is even. But since all the Cfc’s are odd 

integers, this condition becomes 


Mp Ml Mci^ “t“ * * * “1“ Mc^_j^ — 2tii (A9a) 

for some nonnegative integer m. Then from (A8a) 

2p = Mq + Mp + d{2ni — Mp + M 2 + M 2 C 1 + • • ■ + M2cm-i) 

= Mq + 2dTli — [d — I)Mp + C?(M2 + M^2 ci + ■ ■ ■ + M^2cm-l) 

so that the sum Mq + M 2 + • • ■ + M2cm-i is also even: 


Mq + M 2 + • • • + M2cm-i — 2n2, n2 ^ 0. (A9b) 

We define Zi and Z 2 as 

Zi = 277-1 — Mp, Z 2 — 2772 — Mq. 

Then 

2p = Mq + Mp + d(^Zi + Z 2 ), (A9c) 

and we re-express (A9a,b,c) as 

Mi+Mc^-\ - Mc^_^ = Zi, (AlOa) 

M 2 -l- M 2 C 1 -|- • • • -l- M2cm-i — ^2, (AlOb) 

Mq -\- Mp = 2p — d{Zi -|- Z 2 ), (AlOc) 


which formally decouples the three sets of indices. The first two equations require that 0 < Zi,Z 2 < 2m. The 
third equation requires that 0 < Zi -|- ^2 < 2m, and that Zi -|- Z 2 has the same parity as Mq + Mp. Equations (A9) 
however have the additional requirement that the pairs {Zi, Mp) and {Z 2 , Mq) individually have the same parity, so we 
must impose this by hand on the solutions to equations (AIO). (If one set has the same parity, then so does the other.) 


For a given {Zi, Z 2 ), the number of solutions to (AIOc) without this restriction is 
/ 2p+l- d{Zi + Z 2 ) \ 

I . 

If Zi -I- Z 2 is odd, then W is even and there are W/2 compositions (Mq, Mp) in which Mp is even, and W/2 composi¬ 
tions in which it is odd. So to count only the compositions in which Zi and Mp have the same parity, we divide W 
by 2: WW/2, for Zi -h Z 2 odd. 


If on the other hand Zi + Z 2 is even, then W is odd, and there are {W ± I)/2 such compositions for Zi even/odd, 
respectively. 


So, to impose the requirement that {Zi, Mp), (or equivalently {Z 2 , Mq)), have the same parity we make the replace¬ 
ment 

The number of solutions to (AIO) with this requirement is then 


- 

2 


2m 2m —Zi 


KiStxS^)Uo<id,>i=J2 E 2 


Zi —0 Z 2 —0 


2p -|- I — di^Zi -f Z 2 ) -f 


(-1)^1 -f (-1)^= 


Zi -1-777— l\ /Z 2 -1-777— 1 


777 —1 


777 —1 
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From the {Zi,Z 2 ) symmetry of the summand we can replace ( (—l)'^^ + (—1)'^^ )/2 by (—1)'^F Evaluating this 
expression we get 


K{St X 5^^)Uodd.>i = 


4p — 1 


- (p - 1) 


4m + l\ 2m \( 3m — Zi 


2m + 1 / Zi=o 


To evaluate the remaining sum we use the following lemma: 


Lemma 6 


^ 1 

fe=o \ m — 1 


X-k\ ( \X|2^ 

{-If = (-1) — ! 


m — 1 


where [ ] denotes the ceiling function. To prove this identity we write it in the form 


(^) E 

fc=o y m — 1 ^ 

2n+l / k 

(B) e { ^ 

k=o \ m — I 


2n — k\ I n 

1 i' -i\k _ / 1 \m-l I 


(- 1 ) =(-ir 


2n + l — fc\ /n + 1 

(-1)'^ = (-ir-' 


These two identities can be proven by using the identity 


(1 - a;)™+i 


9=0 \ m 


and expressing the product 

(1 +x)™ (1 - 

first as the product of two sums, 


(1 + x)™ (1 - x)'"+i 

then as a single sum, 


\ k=o \ m — 1 


(-x)™-l X™ (-l)'"-lx2™-l(l + x) 

(1 +x)™ (1 - x)'"+i ~ (1 - x2)™+i 


oo / g \ 

E 

g=o \ m / 


= (-ir-'E 


p=o \ m 


(x2P-l +x2P), 


and then equating the even-exponent and the odd-exponent terms on the right-hand side of (All) to those on the 
right-hand side of (A12). 


Now setting k = Zi + m — 1 and X = 4m — 1 we have from this lemma 
2 m / Zi + m — 1 \ ( 3m — Zi\ / 2m \ 

Zi=o \ m — 1 J \ m J \ m j 


And so 


iF(A+xA^^)|dodd,>i = 


4p — 1 


4m \ / 4m + 1 

-{p-f 

2m j \ 2m -I- 1 


Putting all these results together, we get eq.(20b). 
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Appendix B 


We have, from the discussion in the proof of Theorem 3, that 


^o--o 1---1 Ai-'-Ap+i ^i<Zi ••• Apqp;Mi,N 

gi,---,gp 


1 I 9i,---,<7p) 




(Bl) 


(fc) 

To determine the relation between this partition function and ^o-.-o 1---1 consider the case where all 

multiplicities Mq, Mi, ..., M^v-i have no common factor. Then, as discussed earlier, all the cyclic-permutation equiv¬ 
alence classes {ctq = ^p+i, CTi, ..., ctm-i} have N elements. The coefficient is then proportional to 

N times a sum of Q over all possible g’s divided, to avoid overcounting, by the factorials of the multiplicities: 


OC 


M 2 ! • • • Af^v—1 ■ 


Q{k-Aiqi - Apqp-,Mi,N 

qi.---,qp 


1 I qi, ■■■,%) 




(The factor of Ma^^i occurs in the numerator because the multiplicity of Ap+i within the indices {^1 ,..., Ap} is 
only M^p+i - 1.) 


Now consider a particular permutation [cr] = [Ap+i, cti, ... ctat-i] corresponding to a particular set {qi ,..., qp} and 
a particular partition (pi,... , pmi )- [o'] contains, in addition to the Ap+i at the 0 position, an additional — 1 ) 

Ap_|_i’s at various qj positions. The equivalence class containing this [cr], i.e., {cr}, also contains (among others) all 
permutations of the form [cr'j = [ctq = cr^^ = Ap+i, cr},..., (Tq = Ap+i,.. obtained by cycling the ^p+i at 

the qj position to the 0 position, (and the Apj.i at the 0 position to somewhere else). These permutations will also 
be counted in the sum over the q's, under a different set of q's and p’s. Since there are [Ma^+i — 1) such g’s, each 
equivalence class will be counted Maj,+i times in the sum above. We must therefore divide N by an additional Ma^^^ 
to get the number of distinct permutations: 


a(^) 

"^0---0 1---lAi-"Ap+i 


N 

M2\ • • • Ma! — A 


Q{k-Aiqi - Apqp-,Mi,N 

qi.---,qp 


1 191, ■ • ■ ,<?p) 




For example, consider the index set [0011255] and the two equivalence classes below that appear in the sum over 
91,92: 


{5001125} : {pi,p2} = {3,4}; { 91 ,( 72 } = {5,6}; 

{5500112} : {pi,p2} = {4,5}; {( 71 , 92 } = {1,6}. 

Each of the permutations [5001125] and [5500112] are in both of these equivalence classes, so we must divide by 2. 

Now let the multiplicities have a common factor d: gcd{MQ, Mi, M 2 , ■ ■ ■, M^-i) = d. Then there will be some 
equivalence classes of the form {Ap+ia; ■ • - y Ap+ix ■ ■ - y ■ ■ ■ Ap+ix • • • p} for which Oj = cr^+jv/d, so that the “phrase” 
Ap+ia; • • • p is repeated d times. These classes have only N/d distinct elements. The sum in (Bl) will then consist of 
two sums, one over the A^-element classes and one over the (A^/(i)-element classes, which we will express in short-hand 
notation as 


Ak) _ ^Ap+i 

M2\---Mn-i\ 


N 


Ma„ 




N/d 




II 


where division by some factor X is needed to correct for overcounting due to the duplication of equivalence classes 
in the sum. The index Ap+i occurs Maj,+i times in the entire index set [Ap+ix • • • p Ap+ix ■ ■ - y ■ ■ ■ Ap^ix • • • p], and 
MAp+i/d times in the phrase Ap^ix- ■ -y; MA^^i/d is thus the factor X by which each distinct equivalence class is 
overcounted: [N/d)/[M a^+i/ d) =N/Ma^+i- 
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As an example of this situation, we consider [001122333333], with N = 12, Mq = Mi = M 2 = 2, M3 = 6, d = 2, 
and the equivalence classes: 


{303312303312} 

{331230331230} 

{312303312303} 


{pi,P 2 } = {4,10}; {qi,q2,q3,q4,q5,q6,q7} = {2,3,5,6,8,9,11}; 
{Pi,P2} = {2,8}; {gi,(72,<?3,<?4,g5,96,g7} = {1,3,4,6,7,9,10}; 
{Pi,P2} = {1,7}; {91,92,93,94,95,96, <17} = {2,3,5,6,8,9,11}. 


Each equivalence class has N/d = 6 distinct elements, and each class is overcounted by a factor of 3 = M^/d. The 
overall correction factor is then 6/3 = N/M 3 . (Although classes 1 and 3 have the same 9 values, their p values are 
different and will therefore be counted separately by the partition function.) 

Since the correction factor is the same as for the first sum, we can use the same multiplying factor for both without 
any further need to consider them separately: 


Appendix C: Proof of Lemma 1 


For a given set {ki, ..., Kp} let 
j3j = Ki-\ -h Kj, 1 <j <p. 

Then ki = j3i and Kj = (dj — fdj-i for j > 2. In the sums over the k’s in the statement of the lemma, the restriction 
that Ki + Kp < m requires that each j3j < m. And, since each Kj has the range [0, m], each pj-i has the range 
[0,/3j]. We therefore have 


E 


^K1(?1+K292H- 




(kI + • • • + Kp)] 


m 0 p 03 02 

~ ^ ^ ^ ^ ^ ^ . 

/3p=0/3j,_i=0 02=0 01=0 


ljj(.02-0l)q2i^0iqi ^ 


and the statement of Lemma 1 is equivalent to 


N—l p 

n 


m—0 s—1 


1 

1 + 


N—l—p m 0p 03 02 

i-yr E E • • • E E • 

m=0 0p=o0p_i=O 02=0 01=0 



1 ^( 02 - 0 l)q 2 ^ 0 iqi ^ 


The proof of Lemma I, by induction on p, is more direct when stated in this form. 


To prove the lemma we Hrst set p equal to 1 on the right side, multiply that expression by (1 + yuj‘^), and expand: 


N-2 


N-2 


N-2 


m—0 /3—0 m—0 /S—0 m—0 /3—0 

The first sum on the right in eq.(Cl) is 

N-2 N-2 m 

E(-2^r + E(-2/rE^''^ 

m=0 m—1 ^—1 

while the 2nd sum, after relabeling (3 and then m, is 

N — l m 

-E(- 2^rE^''^- 

m=l 0=1 

The 2nd sum in expression (C2) will therefore cancel with all but the m = N — 1 term in (C3), and we have 

N-2 m N-2 N-1 N-1 

m—0 /3—0 m—0 /S—1 m—0 


(Cl) 


(C2) 


(C3) 


(C4) 
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where we’ve used eq.(7) to evaluate the sum over f3 as equal to -1. The lemma is therefore valid for p = 1. 

We now assume the lemma is valid for the product of (p — 1) factors, ns=i(l + and express the product 

of p factors as the difference of two {p — 1 ) products: 


E(-rf"nr 

m—O s—1 






-—In— 

11 /, I -L X 1 _L 7// 


1 + yuj^^ 11 + J x x i y^qs 


ujq^ - ojq^ 


N—p ( 771 04 02 

E ■ E E‘ 

m=l l/9p=0 /32=0/3 i=0 


XPp-Pp-l)qp . . . ,XPi-P 2 )qi, XP 2 -Pl)q 2 , fiiqi 


m fHi fUs I 

-E-EE ^{Pp-Pp-i)qp ,.. ^jj(Pi-03)qi^{P3-0i)q3^j0iqi I ^ 

/3p=0 /33=0/3 i=0 J 

In the first term inside the brackets we make the replacement j32 —>■ Ps, and factor out the gp,..., 74 sums, so we have 


w® - 


m Pa 02 

E-EE ^{0p-0p-l)qp . . .ijj{0A-02)qAi^{02-0l)q2^0iqi 


/3p=0 02=0 01=0 


04 . 03 


-E-EE‘ 

/3p=0 02=0 01=0 


X0p-0p-i)qp ... ,,X04-03)q4, X03-0i)q3,.0iqi 


= V ... V uj^Pp-Pp-^)‘lP . ..^i04-03)q4 ^ a;(fe 0l)q2_^^{03 0l)q3 

a;«2 - w® 

/3p=0 03=0 01=0 


The validity of the lemma for p factors then follows from the identity 


= E 


and a relabeling of the summation indices. A 


Appendix D: Proof of Lemma 2 

To remind ourselves, for a partition Z = {zi, Z 2 , ■ ■ ■, Zj) of p, the vector is defined as 

Qz = 91,92,..., 72, . ,7j,...,7j). 


and the statement of the lemma is that 


E h{(l) = ^ ^ iPlki,k2,...,kp)* h{q^) 


l<qi<---<qp<M 


l<gi T--iqj <M 


where /i is a completely symmetric function of its arguments. 

This proof is also by induction. We separate off the j = p term in the sum and write the lemma in the form 


E k{qi,...,qp) = p! Y ^(9i,...,9p) + E E ^ (p;/ci, fca,..., fcp)* Y (^1) 

'i,...,9p qi<---<qp j=i ki,...,kp qi,---,qj 
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with j = ki + ■ ■ ■ + kp, and the summation limits 1 < Qn M and 0 < < p for all Qn and understood. 

For p = 1 the equality is trivial. For p = 2, this expression reduces to 

h{qi,q2) =2 ^ h{qi,q2) + '^h{qi,qi), (D2) 

91.92 91 <92 91 

which is clearly valid. We now assume that the lemma is valid for summation over the {p— l)-dimensional lattice and 
show that it is then valid for the p-dimensional one. For clarity in the notation, we will denote a partition of (p — 1) 
by F = (t/i, 7/2, ■ • ■, Vj), with multiplicities [/i,..., Ip-i, Ip = 0]. We then have for the summation over p dimensions, 


9l....,9p 


h{qi,.. .,qp-i,qp) 

Qp qi,...,qp-i 

X! Hqi,-■ ■ ,Qp-i,Qp) 


(D3) 


qp 


qi<-"<qp-i 


p-2 


+ Y1 (p- h{ci^,qp) 

j=l Zi....,ip_i 91.•••.93 

by the induction assumption. We now separate the sum over qp applied to the first term inside the brackets into sums 
over the p regions 

(qp < qi < q2 < ■ ■ ■ < <?p-i) (qi < qp < q2 < ■' ■ < qp-i) ■■■ {qi < q2 < ■ ■ ■ < qp-i < qp) 

and the (p — 1) boundaries 

(<7p = <?i < 92 < • • • < 9p-i) • • • (91 < 92 < • • • < 9p-i = qp)- 

By the symmetry of the function h, the p-dimensional sum over each subregion is the same, and we can write (making 
the change p ^ j + \,qp ^ (/j+i in the last sum), 

Y ^(9ii---i9p) = p! Y ^(9i,---,9p) (D4) 

91 .•••.9p 9i<---<9p_i<9p 

+(P-1)! ^ {h{qi,.. .,qp_i,qi) + h{qi,.. .,qp_i,q 2 ) H-+ h{qi,... ,qp-i,qp-i)} 

9l<---<9p-l 

p-2 

+ Y E (-if+Mp-i;^i,^2,...,/p-i,o)* ^ Mqvi99+i)- 

3 = 1 Zi,...,ip_i 9i,...,9j,93 + i 

Now although the functions h{qi,q 2 ,..., qp-i,qi ),... h{qi,q 2 ,..., qp-i,qp-i) in the second term are not completely 
symmetric functions on the (p — l)-dimensional lattice, the sum of these functions is symmetric, and we can apply 
the lemma in its (p — 1) form to it: 

(p-1)! E {h{qi,.. .,qp_i,qi) + h{qi,.. .,qp-i,q 2 ) H-+ h{qi,... ,qp_i,qp_i)} 

9l<---<9p_l 

p-1 

= E E (p-1; ^1,^2,-/p-1,0)* X! {^^x’^(i))+^(qvi9(2))H-h/i(qy,(/(p_i)) }. 

j = l ii,...,ip_i 91.•••,93 

The notation (//j.) here means, for a partition Y, q(^p^ is equal to the value of the fc-th coordinate in the associated 
lattice point . E.g., for Y = (3, 2, 2) and for fc = 2, we have 


/i(9i, 91,91,92,92, 93, 93, 9(2)) = /i(9i, 9i, 9i, 92,92, 93, 93,9i)- 

t t 
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Therefore, 

^ h{qi,...,qp) = p\ ^ h{qi,...,qp) 

qi,...,qp qi<---<qp 

p-1 

+ E E (-If(D5) 

3 — ^ l\ 1 ••• — 1 

X E + hi^Y:<l{2)) + -f ^(qr.'7(p-l)) } 

p-2 

+ E E (-if (p- h{q^,qj+i). 

j=l qi,...,qj+i 

The term on the left and the 1st term on the right are in the “correct” form for the proposition. The 2nd and 
3rd terms on the right contain sums over the partitions of {p — 1) that, to complete the proof, we want to express 
as sums over the partitions of p. Given a partition {yi,... ,yj) of {p — 1), there are two ways a partition of p can 
be generated from it: one can add 1 to an individual part, ym —>■ {ym + 1)) or one can add a 1 to the total sum 
yi + • ■ • + yj ^ yi + • ■ • + yj + ^- As shown below, the 2nd term converts to a sum over the partitions of p by this 
first process, and the 3rd term by the second process. 


We consider the 2nd term first: 

p—1 

E E (-iF”^"^Mp- i;^i,--->^p-i.o)* E E <?(«))• 

i=l ii,...,/p_i qi,...,qj s=l 

For a partition Y = {yi,... ,yj) and for a particular (/(g) such that 

2/1 H-1- Vd-i < s <yi-\ -h 2/d-i + Vd (D6) 

for some d, the lattice point is equivalent to the partition Z = = (yi,... ,yd + 1,... ,yj) of p, since 

/i(qy, <7(5)) = h{qi,... ,qi,q 2 ,... ,q 2 , ■ ■ ■ ,qj, ■ ■ ■, qj,qd), and, by the symmetry of h, we can move the qd in the final 
position to a position among the other qd’s, so there are (yd + 1) qd’s. In terms of multiplicities, if 

Y = [ll, . . . , ly^, ly^J^l, . . . , Ip—I, 0], 

then 


^Vd — [^li • ■ • ) ^Vd 1; ^Vd + l T Ij • ■ • j ^p-l) 0]- 


i / ya I ' j 

In the sum over s, the number of s’s that satisfy condition (D6) above is Pdly^j so the partition Zy^ will appear ydly^, 
times in this sum, and we have, 

p-i 


p—i P~-L 

E E (-i)^”^^Mp- i;^i>--->^p-i,o)* E E 

Zi,...,ip_i qi,...,qj s=i 

p—1 

= E E (-If+'(p-i;/i,---,f ,or^ n;„ ^ Mq. 

j=l h,...,lp-i n=l qi,...,qj 


(D7) 


where the partition has multiplicities [ki,..., kp] given by k^ = — ^m+i = Im+i+^n.m for m = 1,... ,p—1. 

We now multiply and divide on the right by the corresponding multinomial coefficient (p; fci,..., kp)’^: 

p—1 P~1 

^ ^ {p-i;h,...,ip-uor Y E f,?(.)) 

7 = 1 U . qi,....qi s=l 

P-1 


Since 


= 'Z E (-ir‘w fet...... 

j=l h,...,lp-l n=l fP,fcl,---,/Cpj„ 

(p- l]li,...,lp-i,0)* _ 1 n^"~^(^„ - 1)! (n + l)^"+i+^(ln+i + 1)! _ 1 (n + l)(?n+i + 1) 

(p;fci,...,fcp)* p (n + p nl„ 


(D8) 


(D9) 
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we have 

p—1 P~1 

i=l ii,...,/p_i gi....,93-5=1 

= E E (-ir^+^' (p;fci,...,fc,): E Mqx„)- (dio) 

j=l ii,...,ip_i n=l ^ 91.■■■,93 

In this equation, for a given partition K, we sum over all partitions Zn - We now interchange the sums and, for a given 
partition Z, we sum over the corresponding Y partitions. This amounts to dropping the subscripts ‘n' and replacing 
the sums over the Vs by sums over fc’s. So, making the replacement n' = n + 1 and then dropping the prime, 

p—1 P~1 

E E E E 

3 = 1 ;i,...,/p_i 91,...,93 5=1 

= E E {p;ki,...,kpy E 

3 = 1 ki,...,kp n =2 ^ 9i,...,93 

= E E E (Dll) 

3 = 1 fei,...,fep ^ 9i,...,93 

where we’ve used fci + 2fc2 + ■ • • + pkp = p in the last line. 


Turning now to the 3rd term, 

p-2 

E E l;^i.---,^P-i: 0 )* E hici^,qj+i), 

3 = 1 /i,...,/p_i 9i,...,93+i 

the lattice point (q^, ^3+1) is equal to the lattice point for the partition Z = {yi,... ,yj,l) = [^i + l, ^2, • • ■, Ip-i, 0] € 
V{p) with j + 1 parts. Setting j' = j + 1, we can express this term as a sum over the partitions of p from j' = 2 to 

p- 1: 


p-2 

E E i;^i,---,^p-i,o)* E ^(qr,93+i) 

3 — — 1 

p-1 




3 '—2 

p-1 


v,gj + l 

-1 (p- l;li,... ,lp-i,0)* 

(p; /i + 1,..., lp_i, O)" 


(p;li + l,...,/p_i,0)* E 


91,...,93/ 


= E E ^ — ip'^ki,...,kp-i,^)* E ^(^z 

j'=2 fci,...,fcp_i ^ 91,...,93/ 


(D12) 


using 


(p-l;?i,l2,...,/p-i)* ^ 1 1^^+H^i + 1)!2'^?2!---(p-1)^^-^^p-i! ^ h + l ^ h 

(p;/i + 1,... ,lp_i,0)* p 2*2^2! ••• (p - l)'p-iZp_i! p p' 

But since fci = 0 if j' = 1, (corresponding to the partition p = p and kp = 1), the summation lower limit can be 
extended to f = 1 and a “sum” over kp added, along with the replacement (p; fci,..., /cp-i, 0)* —>■ (p; fci,..., kp-i, kp)*. 

Putting (Dll) and (D12) into (D5), the ki/p terms cancel and we get 

p-1 

E Hqi,---,<ip) = p'- E ^(9i>--->'?p) + E E (-1)^"^^”^ (p;^i>--->*p)* E ^(^z)- ^ 

9 i,...,9p 9l<'"<9p i=l fei,...,fep 91.■■■,93 
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Appendix E: Proof of Lemma 3 

The multinomial coefficient (p; fci,..., kp)* has the generating function Q 


J.p 

X! ^ X! {p\ki,k 2 ,...,kp)*x\^ 

V—3 ^15 • • •fkp 


^^2 . . . 


= i:?*‘ 


(El) 


where the sums over the k's are subject to the two constraints 

*1 + ^2 H- 'rkp = j, 

ki + 2k2 + • • • + pkp = p. 

Dividing both sides by jl, summing over j, and interchanging the j and the p summations on the left, this becomes 

p 


tP 


Y. ^ Y {p-,ki,k2,...,kp)* 


p\ 

p=0 j=0 ki,...,kp 

Now setting all Xi to 1, we get 

oo p 


ki ko kn 

•^1 *^2 


1/ ^ \ ry 

= Y.t, E?'* 


j=0 \2=1 


H ^ II Y ^ 1 ’ * 2 ,..., kp)* = exp ^ - = exp (- ln(l - t)) = ^ ^ t". 

p=0 P' j=0 ki,...,kp 

Comparing the two, 


\i=l 


n—0 


Y Y {Plki,k2,... ,kp)* = p\. 

j=0 ki,...,kp 

Now, for a set on nonnegative integers, {/3i, ■ ■ ■, Pq | A > 0}) we apply the differential operator 

001-\-02~\ -|-/3q 

9^1 Xl ■ • • d^lXq 

to both sides of (E2) and again set all Xi equal to 1. On the right we get 


g0i+02-\ —i-/3g 


exp 


9^1 Xl • • • d^lXq 

^ \ t X , 

with no = /3i + 2/32 H-f 9/3g, 

and on the left, 

Q0i+/32-\ -l-/3g IP P 


./3i+2/32 + ---+9/3, “ 1 “ 

^ l/3i . . . g/3, ^ ’ 


l/3i . . . qfi 


n—0 


Y ~Y Y {P',ki,k2,...,kp)*x\^ 

p=0 P' j=0 ki....,kp 


^k2 ^ . ^kp 


Xi — 1 


d^^Xl- ■ ■ d^lXq 2-^ p\ 

P — 0 

^ fp 

= Y ~\Y Y ■ ■,fcp)* fci(fci -1) ■ ■ • (fci -/3i +1) 


oo P 


pi 

p=no j=0 ki,...,kp 


oo p 


X k2ik2 - 1) • • • (/i:2 - /32 + 1) ■ • ■ kq{kq - 1) • • • {kq - Pq 1) 
kn 


= Y hY Y 

p=no ^ j=0 ki,...,kp 


ki 

Pi 


■•/3q!( I {p;ki,k2,...,kp)*. 

Pq 


Therefore 


P [ ki 

^ ^ L 

j=0 ki,...,kp y Pi 

for p > /3i + 2/32 H-1- pPp- A 


{p]ki,k2,...,kp)* = 


p\ 


1A/3i! 2/32/32! •••p/3p/3p! 


(E2) 


(E3) 


(E4) 


(E5) 


(E6) 


(E7) 
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Appendix F: Sample calculations 


As examples, we will use Theorem 3 to calculate the N = 7 coefficient (70123456 and the N = 8 coefficient (7oo224466- 


1- (7 o 123456 


Since the Ak’s are all distinct in this example, the set of partitions of the set {Ai, A2, A3, A4}, i.e., 'P{AiA2A3A4}, 
and the set of partitions of the set of indices of the A^’s, 7^{1234}, have the same number of elements. We will 
therefore use P{1234} and the somewhat simpler expression, (lOd), to find (7oi23456: 


(7o123456 


T’{1234} 


-7x<'5!+ ^ ^ H 

ee'P{1234} k=l Ai....,Aj=0,l 

5 - ELi MX{0t) + Zt) \ J ( X{9s) + Zs-1^ 




1 - ELi ExiOt) 


n 


2:.= — 1 


{ (1234), (123)(4), (124)(3), (134)(2), (234)(1), (12)(34), (13)(24), (14)(23), 
(12)(3)(4), (13)(2)(4), (14)(2)(3), (23)(1)(4), (24)(1)(3), (34)(1)(2), (1)(2)(3)(4)} . 


The value of, for example, X{9 = 1234) is Ar(1234) = — (Ai + A2 + A3 + A4) mod 7 = —(2 + 3 + 4 + 5) = 0. The rest 
of the X{9)'s are: 


A:(123) = 5, A:(124) = 4, X(134) = 3, A(234) = 2, 

A(12) = 2, A(34) = 5, A:(13) = 1, A:(24) = 6, A:(14) = 0, A:(23) = 0, 

A:(1) = 5, X(2)=4, A(3)=3, A:(4) = 2. 


The subset of 7^{1234} of partitions that correspond to non-zero terms, for which 1 — X) EX{9t) > 0 for at least some 
values of the At’s, is then: 

{(1234), (13)(24), (14)(23), (13)(2)(4), (14)(2)(3), (23)(1)(4)} . 
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So 


Coi23456 = -7 { 5! + (-7)^^3! 


5 -4Ai 


Ai = l 


e=(1234) 




5 - 3Ai - 8 A 2 


1 — Ai — 6 A 2 




+ (-7)^i+^n!l! 


+(-7)^i+^"+^n!0!0! 


_|_ ( —7)'^l+'^2 + '^3 J^IQIQI 


_|_(_7)'^1 + '^2+A3]^|Q|Q| 


1 — oAi — 0 A 2 y Y 1 
5 - 3Ai - 5 A 2 - 3 A 3 


1 — Ai — 4 A 2 — 2 A 3 

5 - 2Ai - 5 A 2 - 4 A 3 
1 — OAi — 4 A 2 — 3 A 3 

5 - 2Ai - 6 A 2 - 3 A 3 
1 - OAi - 5 A 2 - 2 A 3 


5 - 2Ai - 2 A 2 \ / 1 \ V 1 


A 2 


A 2 


(Ai.A 2 ) = (1,0) 


e=(13)(24) 

(Ai.A2) = (1,0),(0,1),(14) 


e=(14)(23) 


2 \ ^ V 4 \ ^ V 2 ^ 


1 \ / 4 \ / 3 ^ 


1 y V 0^ 

1 y ^ / 5 \ ^ V 2 ^ 


(Ai,A 2 ,A 3 ) — ( 1 , 0 , 0 ) 


e=(13)(2)(4) 
(Ai,A 2 , A 3 ) —(1,0,0) 


e=(14)(2)(3) 
(Ai.A2.A3)) = (1,0) 


e=(23)(l)(4) 


where for clarity we’ve indicated the particular partition that each term corresponds to and the allowed values of the 
At’s. We have then 


C'oi23456 — — 



+ (-7)" 







-105. 


2. (700224466 


In this example Mi = 0 and the non-zero terms are characterized by 'Y^\tX{6t) = 0. So only partitions 0 = 
(di) • • • {9j) that contain at least one “null” part 9k such that X{9k) = 0 , will contribute; for a non-null part 9q, Xg 
must equal zero for the term to contribute to the coefficient. When Mi = 0, the two binomial coefficients in the 
formula, 


N-Mo-l-EMXm + zt) 
Mi-J2XtX{9t) 


and 
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can both be replaced by 1. The first binomial coefficient equals 1 since, for non-zero terms, it has zero in it’s lower 
entry; in the second coefficient, either X{9s) = 0 or As = 0. So we have 


Co0224466 — 


5! 


2 ! 2 ! 2 ! 


ee'P{22446} 8 ” k=l 


{zk - 1)! 




X ^(-8)^ 


E 


-hAj H 


-^XtXiOt 


The null subsets of {2,2,4,4,6} are (2446), (224), (44), and (26). The 10 contributing partitions are then: 

(2446)(2): (Ai, A 2 ) = (1,0) 

(224)(46): (Ai, A 2 ) = (1,0) 

(226)(44): (Ai, A 2 ) = (0,1) 

(244)(26): (Ai, A 2 ) = (0,1) 

(224)(4)(6): (Ai, A 2 , A 3 ) = (1, 0, 0) 

(44)(26)(2) : (Ai, A 2 , A 3 ) = (1, 0 , 0 ), ( 0 , 1 , 0 ), ( 1 , 1 , 0 ) 

(44)(22)(6): (Ai, A 2 , A 3 ) = (1, 0, 0) 

(26)(24)(4): (Ai, A 2 , A 3 ) = (1, 0,0) 

(44)(2)(2)(6) : (Ai, A 2 , A 3 , A 4 ) = (1, 0 , 0 , 0 ) 

(26)(2)(4)(4) : (Ai, A 2 , A 3 , A 4 ) = (1, 0 , 0 , 0 ) 


(To simplify the expression below, we will omit the factor ]([g when all the k’s equal 0 or 1.) Summing over these 
partitions, we get 

(1^00224466 = —i 


15- 


(4-l)!(l-l)! (3-1)! (2-1)! (3-1)! (2-1)! (3 - 1)! (2 - 1)! 


1 ! 2 ! 1 ! 1 ! 0 ! 0 ! 


4-2 X 


( 2 - 1 )! ( 2 - 1 )!( 1 - 1 )! 


2 ! 1 ! 0 ! 0 ! 1 ! 1 ! 2 ! 0 ! 1 ! 0 ! 2 ! 0 ! 
(3-1)! (1-1)! (1-1)! 

2 ! 1 ! 0 ! 0 ! 1!0 0 ! 0 ! 1 ! 

(2-1)! (2-l)!(l-l)! 


1!2!0!! 1!0!1! 


(2-1)! (2-l)!(l-l)! 


0 ! 2 ! 0 ! 1 ! 0 ! 1 ! 1 ! 0 ! 0 ! 0 ! 2 ! 0 ! 2 ! 0 ! 0 ! 0 ! 0 ! 1 ! 1 ! 0 ! 1 ! 1 ! 1 ! 0 ! 0 ! 1 ! 0 ! 

1 (2 - 1 )! (1 - 1 )! (1 - 1 )! (1 - 1 )! 1 (2 - 1 )! (1 - 1 )! (1 - 1 )! (1 - 1 )! 


1!2!1! 0!2!0! 1!0!0! 1!0!0! 0!0!1! l!l!2! 1!0!1! 1!0!0! 0!1!0! 0!1!0! 

64 (2 - 1)!(2 - 1)!(1 - 1)! 

T 

32 


0 ! 2 ! 0 ! 1 ! 0 ! 1 ! 1 ! 0 ! 0 ! 

= -120 -b 96 -b 32 -4 16 -f 32 -b 32 -b 32 -4 8 -f 32 -b 8 -b 16 - 128 = 56. 


For this coefficient, the W-J method is actually more efficient. The null subsets of {2,2,4,4,6,6} are (224466), 
(2266), (2446), (466), (224), (44), and (26). There are then 5 null multisets to be summed over: (224466), (2266) (44), 
(2446) (26), (224) (466), (26) (26) (44). We have then, using Wyn-jones’ formula, 

= (-8)1^ . (-8)- . (-84 + (-8)4^ - 


2 ! 2 ! 2 ! ' ' 2 ! 2 ! 2 ! 
3^ (2-l)!(2-l)!(2-l)! 
2 ! 2 ! 

= -120 -b 48 -b 192 -b 64 - 128 = 56. 


2 ! 


2 ! 2 ! 


Appendix G: The coefficients for the N — 6,7,8 determinants 

In this section we give our results for the coefficients [ {C'mo'"Mn_i } ] N=Q, 7, 8, using Theorem 3. The results 
are expressed in the form of a table. In the first column are listed the allowed partitions of N; the second and third 
columns list, respectively, the form of the coefficient Cao---aN-i and all the additive multiplets associated with this 
partition; and the forth column lists the values of [ {C'mo...mn_iI ]■ 
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The subscript n on is used to indicate the number of elements in a particular additive multiplet. 

To simplify the notation however it is understood that if there is no subscript then the multiplet has N elements. 

Our results are: 


N = 6: 

(12 super 

-multiplets). [ ] 


6 

C'oooooo 

{C'eoooool 

1 


42 

^OOOOaa 

{^ 400200 } 

-3 


411 

C'oOOOab j 

f {C'lioooi} J 

[ {C'loioiol 1 

f -6 

1 -6 


33 

C^OOOaaa 

{C' 303000 } 

2 


321 

CoQOaab 

{^ 320010 }’ {C' 301002 } 

6 


3111 

CoQOabc 

{C' 311100 }) {C' 300111 } 

12 


222 

CoQaabb j 

f {C' 222000 } 

L {C'20202o}n=2 

f -9 

1 9 


2211 

CoQaabc j 

^ {C' 210120 } J 

L {C' 210201 } 1 

f -18 

1 0 


21111 

CoQabcd 

{C' 211011 } 

0 


N = 7 : 

(12 super 

— multiplets) 



7 

C'ooooooo 

{C'tooooooI 


1 

511 

C'oooooab 

{^510000l}j {*^501001o}> {C' 5001100 } 

- 7 

421 

^OOOOaafc 

{C'420001o}> {C 4021000 }) {C'4io200o}j {^*4000201 

}> {^’ 4000120 }’ {C' 4010002 } 

4111 

CooQQabc 

{C'4iioiOo}> {^ 4001011 } 


14 

331 

CoQQaaab 

{C' 3300100 }) {C'sisooooIj {C' 3003010 } 

- 7 

322 

CooQaabb 

{C'3200002}> {C' 3020020 }) {C' 3002200 } 

14 

3211 

CooQaabc 

1^3211000/5 1^3020101/5 1^3012001/5 1^3100210 

}) {C' 3101020 }) {C' 3000112 } ~ 21 

31111 

CooQabcd 

{C' 3110011 }) {C'sohhoIj {C' 3101101 } 

7 


2221 

22111 

1111111 


C^Qaahbc 

C^Qaabcd 

C^ahcdef 


{^1002022}) {*^1220200} 

{C'roi221o}> {^'1201102}) {1^1120021} 
{C'ri0220l}> {C'r210012}j {1^1021120} 

{C'lllllll}n=l 


- 7 

-14 

35 

- 105 
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N = 8: 
8 

62 

611 

521 

5111 

44 

431 

422 

4211 

41111 

332 

3311 

3221 

32111 

311111 

2222 

22211 

221111 

2111111 


(49 super 
Coooooooo 

^OOOOOOaa 

C'oOOOOOah 

^OOOOOaab 

^OOOOOabc 

^OOOOaaaa 

^OOOOaaab 


Q 


OOOOaabb 




OOOOaabc 


GoOOOabcd 


Q 


OOOaaabb 


GQQQaaabc 


G 


OOOaabbc 


GQQQaabcd 

Goooabcde 

C^Qaabbcc 

GoOaabbcd 

G^Qaabcde 

G^Qabcdef 


multplets) 

{^80000000 } 

{^60002000} 

f {C'eioooooilj {C'eooioiool 
1 {C'eoioooio} 

I {^ 52000010 }> {^ 50120000 }> {^ 50000210 }> {^ 50100002 } 
\ {C' 5020100 o}> {C'50001020} 

/ {^ 5110010 o}> {^50010011} 

1 {C'sioiioooli {C'50001101} 

I {C4040OOO0} 

\ {^40004000 }"=4 

{*^43000100}) {C' 4003000 l}> {^ 4100030 o}> {C4001OOO3} 

{C 4202000 o}> {C40000202} 

{C40200020} 

{C4200OOO2}) {C40020200} 

{C42101000}> {C4oo2101o}> {C4oiol20o}> {C4000IOI2} 
{C 41210000 }> {C 4 ioio 02 o}> {C 4020010 i}> {C4000OI21} 
{^41002001 }> {^ 40012100 } 

{C40102010} 

{^ 4110001 l}j {*^40110110} 

{C'41010101} 

{C'IiooiiioIj {C'40111001} 

{ C'23000003 } J { C'20030300 } 

{C'20300030} 

{C'33110000}) {C'31030010}) {C'30100301}) {C'30000113} 
{*^00130310} 

{*^01030301} 

{*^01300031} 

{*^3210002o}> {^3022001o}> {C'3010022o}) {C'30200012} 

{C'3200120o}> {C'30021002} ] 

{C' 3200201 o}> {C' 3 oi 2200 o}> {C' 3 oo 0221 o}> {C'30102002} 

{C' 32 l 0010 l}> {C' 3002011 i}j {C'31110200}) {C'31010012} 
{C'3201011o}> {(^31120001}, {(^31000211}, {(^30110102} 

{^ 3201100 l}> {^ 3102110 o}j {^ 3001120 l}j {^ 31001102 } 
{C' 3 i 20100 l}> {C' 3 oo 1112 o}j {C'30211100}j {C' 3100102 l} 

{C'3120011o}> {(^30110021} 

{C' 3 iio 210 o}> {C'30012011} 

{C'siiiioiolj {C'30101111} 

{*^20022200} 

{*^02200022} 

{ *^20202020 }n =2 

{C'ooi2221o}> {C'02102012} 

{C'oi 02220 l }5 {C'02012102} 

{^ 2212001 o}> {^20100212} 

{C'oo21212o}> {C'01202021} 

{^2211110o}> {^2102101l}> {^2110120l}> {^20011112} 
{C'02110112} 

{C'21102011} 

{C' 2101210 l}n =4 

{C'21110111} 


[ {Gldg-.-Mj} ] 
1 

-4 


16 

16 

-2 

6 


-12 

20 

20 


-24 

-24 

8 


8 


16 

16 

-48 


16 

-16 


32 

32 

-32 

-32 

-16 

48 

16 


-32 

32 

-32 

32 

-32 

32 

64 

8 

56 

56 


{ -48 
16 
-80 
-16 

{ -32 
32 

-160 

96 

-64 
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To avoid writing out all N of the circular permutations for the N = 7 and 8 determinants, we will use the short-hand 

notations “ H- " and “ ± • • •" to mean, respectively, “plus all cyclic permutations of the preceding terms” and “plus 

all alternating-in-sign cyclic permutations of the preceding terms”. Then the determinants are: 


dei[A,B,C,D,E,F] = 

- 3 [ - B^E^ + + E^B"^ - ] 

- 6 [ A'^BF - B'^CA + C'^DB - D^EC + E^FD - F^AE ] 

- 6 [ A'^CE - B^DF + C'^EA - D^FB + E^AC - F^BD ] 

+ 2 [ A^C^ - B^D^ + C^E^ - D^F^ + E^A^ - F^B^ ] 

-h 6 [ A^B'^E - B^C^F + C^D^A - D^E'^B -h E^F'^C - F^A^D ] 

+ 6 [ A^CF"^ - B^DA^ + C^EB^ - D^FC'^ + E^AD^ - F^BE"^ ] 

+ 12 [ A^BCD - B^CDE + C^DEF - D^EFA + E^FAB - F^ABC ] 

-h 12 [ A^DEF - B^EFA + C^FAB - D^ABC + E^BCD - F^CDE ] 

- 9 [ A^B'^C'^ - B'^C'^D'^ + C'^D'^E'^ - D'^E^F^ + E'^F^A^ - F'^A^B'^ ] 

+ 9 [ A^C^E^ - B^D^F^ ] 

- 18 [ A^BDE'^ - B^CEF^ + C^DFA^ - D^EAB^ + E^FBC^ - F^ACD^ ]; 

det[A,B,C,D,E,F,G] = A’^ + B'^ + C’^ + D’^ + E’^ + F’^ + C’^ 

- 7 [ A^BG + A^CF + A^DF -h • • • ] 

-h 7 [ A'^B^F + A^C^D + A^BD"^ + A*E^G + A'^FF^ + A^CG^ -f • • • ] 

-h 14 [ A'^BCF + A^DFG -f • • • ] 

- 7 [ A^B^F -h A^BC^ + A^D^F -h • • • ] 

-h 14 [ A^B'^G'^ + A^C'^F"^ + A^D'^F'^ -h • • • ] 

- 21 [ A^B^CD + A^C^FG + A^CD'^G + A^BE^F + A^BDF'^ + A^FFG^ -f • • • ] 
-h 7 [ A^BCFG F A^CDEF + A^BDFG -f • • • ] 

- 7 [ A^B^CF^ + A^B^D^G -f • • • ] 

- 14 [ A^B^CEG + A^BC^FF + A^D^FFG -f • • • ] 

-h 35 [ A^B^DFF + A^BC^DG + A^BCD^F -h • • • ] 

-105 ABCDFFG- 
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det[A, B, C, D, E, F, G, H] = A® - + C® - Z?® + S® - F® + G® - iZ® 

- 4 [ ± •••] 

- 8 [ A^BH + A^DF ± • • • ] - 8 [ A^CG ± ■ • • ] 

+ 8 [ A^B'^G + A^GD^ + A^F^G + A^GH^ ± • • • ] + 8 [ A^G^E + A^EG^ ± • • • ] 

+ 16 [ A^BCF + A^DGH ± • • • ] + 16 [ A^BDE + A^EFH ± • • • ] 

- 2 [ A^G"* ± • • • ] + 6 [ A'^E'^ - B^F^ + C^G^ - D'^H^ ] 

- 8 [ A-^B^F + A^D^H + A^BF^ + A^DH^ ± • • • ] 

- 12 [A^B'^D'^ + A^F^H"^ ±---] + 2Q [ A^G^G^ ± • • •] + 20 [ A^B^H"^ + A^D^F^ ± ■ ■ ■] 

- 24 [ A'^B'^GE + A^D'^EG + A^CEF"^ + A'^FGH^ ± • • • ] 

- 24 [ A^BG^D + A'^BDG^ + A^C'^FG + A^FG^H ± • • • ] 

+ 8 [ A^BF^H + A^DE'^F ± • • • ] + 8 [ A^GE'^G ± • • • ] 

+ 16 [ A^BCGH + A'^GDFG ± • • •] + 16 [ A'^BDFH ± • • •] 

- 48 [ A'^BEFG + A'^GDFH ± • • • ] 

+ 16 [ A^B'^G^ + A^D^G^ ± • • • ] - 16 [ A^G'^E^ ± • • • ] 

+ 32 [ A^B^GD + A^BD^G + A^CF^H + A^FGH^ ± • • • ] 

+ 32 [ A^C^DH ± • • • ] - 32 [ A^G^FG ± • • • ] - 32 [ A^BDF^ ± • • • ] 

- 16 [ A^B'^GG^ + A^C^D^G + A^GF'^G'^ + A^C^GH^ ± • • • ] 

+ 48 [ A^B^EF"^ + A^D^FH^ ± • • • ] 

+ 16 [ A^B^F^G + A^CD^F^ + A^F^F^G + A^CF^H^ ± • • • ] 

- 32 [ A^B^GFH + A^D'^FGH + A^BCDF^ + A^BDGH^ ± • • • ] 

+ 32 [ A^B^DFG + A^BCD'^H + A^BF'^GH + A^CDFH^ ± • • • ] 

- 32 [ A^B^DFH + A^BD^EF + A^DEF'^H + A^BEFH'^ ± • • • ] 

+ 32 [ A^BG'^EH + A^DEFG'^ + A^C^DFF + A^BFG^H ± • • • ] 

- 32 [ A^BG'^FG + A^GDG^H ± • • • ] 

+ 32 [ A^BGE'^F + A^DF^GH ± • • • ] 

+ 64 [ A^BCDFG + A^CFFGH ± • • • ] 

+ 8 [ A^B^G'^F'^ ± • • • ] + 56 [ A^B'^D'^E'^ ± • • • ] + 56 [ A^G'^E'^G^ - B^D'^F'^H'^ ] 

- 48 [ A^B^G'^DH + A^BD^FG^ ± • • • ] + 16 [ A^B'^C'^EG + A^GD^EG"^ ± • • • ] 

- 80 [ A^B^GD^G + A^CF'^GH'^ ± • • • ] - 16 [ A^BG'^DE^ + A^BDE'^G'^ ± • • • ] 

- 32 [ A^B'^GDEF + A^BD^FGH + A^BCFF^H + A^DFFGH^ ± • • • ] 

+ 32 [ A^BGFFG^ ± • • • ] - 160 [ A^BCE^GH ± • • • ] 

+ 96 [ A^BDE^FH - B^CFF^GA + G^DFG'^HB - D^FGH^AG ] 

- 64 [A^BCDFGHF---]. 

There are various methods that can be used to check these results. One way is to use the identities 

det[l,l,...,l]=0, 
det[0,l,...,l] = (-l)'^"'(iV-l). 

[The second identity follows directly from eqs.(2) and (7)]. The first identity does not however constitute a check 
for even TV; since the coefficients in additive multiplets alternate in sign, it is satisfied automatically. As sums over 
multiplets, these identities are: 


{C'moMi-.-Mjv-i} ] “ (TV odd), 

N-l 

fc =0 
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Another method, when N is divisible by some integer d and when the coefficients for the smaller (A^/d)-dimensional 
determinant are known, is to use the relation 

det[a;o, 0,..., 0, 0,..., 0, X 2 d, 0,..., 0, XN-d, 0,..., 0] = ( det[xo, Xd, X 2 d, • ■ ■, XN-d] )‘^, 

where the only non-zero matrix elements in the determinant on the left are those of the form Xnd mod n- For example, 
for = 6 and d = 2, 

det[A, 0, C, 0, A, 0] = ( det[A, C, E] f = {A^ + - iACEf . 

Expanding the right side, all coefficients of the form C'moOM 20 M 40 ^le read off and compared with the calculated 
values. 
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